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PREFACE 






This report describes a multi-state project developed and coordinated 
by the Federation of Rocky Mountain States. The project developed and 
tested methods for combining earth resources satellite (lANDSAT) data 
with other multi-source data via computer mapping techniques for use in 
natural resources and land use planning in the Western States - 

The project was carried out by mutually defined procedures in six states — 
Montana, Wyoming, Colorado, New Mexico, Utah, and Arizona. Colorado State 
University and Los Alamos Scientific Laboratory provided technical assistance. 
Two interstate areas of 5,000 square miles each and two intrastate areas of 
3,000 square each were delineated for all subsequent LANDSAT and high altitude 
remote sensing during the project. Four 7.5-minute quadrangles x^^ere selected 
wi thing these large areas as test sites. LANDSAT computer-compatible tape 
classification mapping for 1.1 acre cells x^as conducted with multiple date 
imagery. Land use and cover categories were selected by the states, ranging 
from 19 categories in one state to 81 in another. 

In order to place the LANDSAT application within the context of a regional 
data bank, various non-LANDSAT maps were collected on complementary topics. 
These were all converted into the 1.1 acre grid system and computer composited 
x^ith the LANDSAT maps for deriving and displaying the complex patterns in- 
volved in determining feasibility for surface mining or urban development, etc. 
A key purpose was to demonstrate the appropriate mix of LANDSAT utilization, 
data banking and compositing relevant to regional planning. 

The approach was designed for large areas of interspersed federal, state, 
and private lands, with dynamic interrelationships in mineral, x^^ater, agri- 
cultural and recreational land uses. This approach raised numerous administra- 
tive questions, such as standardization of the ground truth data analysis, 
standardization of land use categories/subcategories, and systematic central 
processing of LANDSAT tapes. 

For this project, state governments designated state lead agencies to 
coordinate work among other state agencies. Through the lead agencies, state 
participation and understanding of LANDSAT use and regional data banking 
progressed toward more centralized operations in most participating states. 

Although this was a relatively brief and modest project in this large and 
dynamic region, it developed several innovative procedures x?orthy of continu- 
ation: (1) a new LANDSAT Mapping System (LMS) for LANDSAT digital interpre- 

tations; (2) use of cellular mapping and compositing for combining the LANDSAT 
maps with other forms of data; (3) training new groups in each state; (4) 
defining the need for a ground truth manual; and (5) the need for fixed, repeti- 
tive ground truth sites for continual LANDSAT use. 

This project marks the beginning of state, interstate, and federal collabor 
tion in these techniques which should now be converted into more extensive opera 
tional systems to meet the characteristic problems of energy, agriculture, 
settlement, and water utilization in the West. 
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1 . 0 INTRODUCTION 


1.1 Scope 

This is the Final Report of a multistate project to demonstrate the 
utility of satellite remote-sensed data in combination with other multi- 
source data as an aid to natural resources and land use planning 
decision-making. 

This was a genuinely cooperative regional project involving the 
six Rocky Mountain states of Arizona, Colorado, Montana, New Mexico, 
Utah, and Wyoming. Colorado State University (CSU) and Los Alamos 
Scientific Laboratories (LASL) provided technical assistance, while 
the Federation and its Earth Resources Technology Applications Com- 
mittee of Its Natural Resources Council coordinated all participants. 

This effort combined LANDSAT digital data with other types of data 
from a variety of sources via a computer composite mapping process 
(see Figure 1.1). The process and results are designed to demonstrate 
the utility and potentials of such multisource analyses for use in 
solving land use and other geographic-related planning and decision- 
making problems. 

1.2 Background 

This project is closely related to and logically follows past 
Federation projects that developed and demonstrated the use of computers 
for multisource mapping* and geographic analyses.^ The U.S. Economic 
Development Administration (EDA) has sponsored continuing Federation 
efforts in this area. 


1t^ Federation of Rocky Mountain States; also, FRMS. 

2 

See Section 7.0 — References/Bibliography — for a complete listing 
of FRMS reports on these subjects. 
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The need for this type of utomated data analysis (and storage) 
system is mounting, especially in the West. Many new and pressing de- 
mands are being made on the region's resources and way of life. De- 
dicionmakers from all levels of government and the private sector are 
faced with increasingly complex problems concerning land use, resource 
management, envxroiuaental quality, and socioeconomic development. 

Issues of energy development, urban planning, air quality, agri- 
culture production, water allocation, and so on, are compounded because 
of inter j urisdictional authority, mixed land oxmership patterns seasonal 
fluctuations, and other common situations. Attention to these issues 
on the part of the states, the federal government, private business and 
regional organizations within the large Western area, has created an 
enormous amount of data and Information that can be used to make 
decisions — at both the policy and program levels. Therefore, a system 
that can analyze this variety of multisource data and show the implica- 
tions of various decisions is highly desirable. 

NASA.^S'^ LANDSAT earth resources satellite is a new, powerful source 
of information that can contribute substantially to more informed decision- 
making. The LAMDSAT is continually circling the earth, returning to the 
same point every 18 days. It is producing multispectral data that can 
provide a user with information at 1 ,w cast for large areas. This informa- 
tion can be in the areas of geology, land cover, land use, water quality, 
crop conditions, etc. 

This project demonstrated the repeatable LANDSAT survey procedures, 
together with composite mapping of all other types of data "invisible" to 
LANDSAT but needed for resource planning. It placed remote sensing and 
particularly LANDSAT digital data into its appropriate role in the overall 
process of producing complex surveys of agriculture, strip mining, urban 
development and water management, where nxmierous other forms of evidence 
in addition to LANDSAT data are needed on land ownership and leasing, 
industrial activity, subsurface resources, etc. 

1.3 Objectives 


The project was aimed at several long-term objectives of the Federa- 
tion of Rocky Mountain States: (a) improved land and water use studies 

through systematic surveys, automated data banking and analysis; and (b) 
the transfer of technologies between levels of government, and especially 
to improve the staters capabilities to make better decision 

Specific objectives were to: 

o encourage interstate cooperation in the utilization of 
earth resources satellite technology for solving land 
use planning problems; 


3 


National Aeronautics and Space Administration. 
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• \vork toward the development o£ compatible inter- 
agency, interstate info rmation system' procedures*} 

• test and adopt a common land use classification; 

o evaluate the efficiency and applicability of using 
satellite data in combination with other data in a 
land use information system; and 

e provide a mechanism for the exchange of information 
about remote sensing and geo-information systems. 

1.4 Approach 

To accomplish these objectives a series of demonstrations in each 
of the six participating states was initiated. A lead agency in each 
state worked with the Federation and -the technical assistance organiza- 
tions following these general tasks: 

« Procurement and^.preparation of data for later LANDSAT data 

analysis, including: LANDSAT computer compatible tapes (CCT’s); 

airplane coverage of the test sites (see Figure 3,1) ..and ground- 
truth data. 

• Identification and verification of land use and land cover 
information via statistical signature analysis of LANDSAT 
digital data, and field verification. 

a Mapping, via computer, of the land use/cover data at a 
1.1 acre resolution (cell size) and a scale of 1:24,000. 

a Collection of raultisource ancillary data for combination 
with LANDSAT data. Selection by each state of appropriate 
land use analysis problems and data values. 

• Computer composite mapping of the LANDSAT land use/cover 
data with the other state-collected data to produce composite 
maps indicating solutions to complex land-use and energy- 
planning problems. 

• Evaluation of the results, communications through meetings 
and reports and final report production and dissemination. 

The Federation coordinated and managed the project and is respon- 
sible for the compilation of this project report, parts of which were 
written by all participants. 

The Federation's Earth Resources Technology Applications Committee 
provided broad-based advisory services to project participants. 
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Each state lead agency handled ground truth, ancillary data 
collection and compositing aspects. 

CSU conducted all LANDSAT data processing, using newly developed 
third- generation software programs; LANDSAT ^fapping System CLMS) . 

LASL provided the bulk of multisource computer compositing. Utah 
and Montana utilized their own computer mapping programs. 

In summary, this project is one of the most extensive geographic, 
interstate applications spanning a variety of test sites from a northern- 
border state to a southern-border state, from forest to desert, and from 
mining to urban development uses, with the participation of many state 
agencies . 

1.5 Summary Recommendations 


Regional Data Analysis Needs : 

« Additional investments are required in designing and 
developing institutional mechanisms for increasing 
earth resources technology utilization. 

® Regional data banking and analytic services should 
include a variety of multisource data and formats, 
be flexible for serving many users, and cover large 
areas rapidly. 

• A regional earth resources technology applications 
center should be established in the West to provide 
for user awareness and education, information exchange 
and networking, technical assistance and capacity 
building, and mobilization of resources within the 
region. 

For Centralized State Services : 

o Additional information, training, and assistance should 
' be given state agencies corresponding with their present 
capabilities and needs for multisource data analysis. 

® Upgrading of LANDSAT products and use with other multi- 
source data will encourage state efforts in central data 
use and managment 

• Standardized procedures for using LANDSAT, improved accuracy 
and low cost will also assist in developing state efforts. 
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For Operational Use of LANDSAT : 

e Improved LANDSAT products such as computer-enhanced 
imagery and composites of LANDSAT with base maps or 
other data will provide users with more acceptable 
results. 

• A standard training site procedures manual for use by 
nontechnical personnel would assist users in efficient, 
easy utilization of LANDSAT by streamlining ground truth 
and verification procedures. 

• A series of standard, permanent training sites in each 
state and designation of a lead state monitoring agency 
would help to take LANDSAT use out of the project stage 
and into the operational stage by supplying immediate 
and detailed ground truth data. 

The LANDSAT Mapping System (LMS) : 

• Optimum use of LANDSAT through automated digital analysis 
procedures requires users understand and supply adequate 
information for class selection, training data selection, 
and signature development. 

• Documentation of the LMS system and its use are required 
for efficient, continued use of the system. 
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2.0 FINDINGS AND RECOMMENDATIONS 


2.1 Regional Data Analysis Needs 


2.1.1 Findings 


In the sparsely settled and scenic states of the Mountain 
West, land use and natural resource issues are especially impor- 
tant. Increased energy development and potential boomtowns, com- 
bined with the region's water, agricultural and recreational needs, 
make these issues recurring topics in to^m councils, state legis- 
latures and governors ' offices throughout the Rocky Mountain States . 

Businessmen, too, are often faced with land-use decisions, 
including where to locate industries and plants for optimum pro- 
ductivity and profitability — with the least disruption to the 
environment and the most benefit to the people nearby. 

Public concern about the use of some land may be clear now, 
but areas of future concern and land-use conflict require more 
than a crystal ball to deduce. Decision-makers need to consider 
all components of the present and future environment in order to 
make responsible plans. 

> 

Some o£ the primary land and resource management programs 
mandated by state and federal laws are:^ 


• Land Use Planning • 

• Wetlands Management • 

• Coastal Zone Management * 

• Floodplain Management • 

• Water Quality Management • 

• Agriculture • 

• Fish and Wildlife • 


Transportation Planning 

Forestry Management 

Water Resources Planning 

Land Reclamation 

Air Quality Management 

Solid Waste Management 

Environmental Impact Statements 


^National Conference of State Legislatures, 1976. Final Report of 
the Task Force on State Use of Satellite Remote Sensing, Denver, CO 
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ORIGINAL PAGE IS 
OP POOR QUALITY 

These activities require a variety of data types, such as 
shown in the following categorization:^ 

CATEGORIES, AND SUBCATEGORIES OF BASIC DATA 


1. Cartographic Data 4. Water Resources 


a. 

Line Maps 

a. 

Surface 

b. 

Low Altitude Photography 

b. 

Subsurface 

c. 

High Altitude Photography 

c. 

Man's Activities 

d. 

Satellite Imagery 



e* 

Special Purpose Data 

5. Geologic and Land Resources 

Meteorological 

a. 

Surface 



b. 

Subsurface 

a* 

Climate 

c* 

Man's Activities 

b. 

Air Quality 



c. 

Man*s Activities 

6. Socio-Economic Resources 

Biological Resources 

a* 

Social 



b. 

Economic 

a* 

Animal 

c* 

Commerce 

b. 

Plant 

d. 

Governmental 

c. 

Ml cro- organisms 

e. 

Archaeological 


Many other studies 3 have demonstrated this need for a variety 
of wide- reaching types and scales of data from a variety of sources. 

NASA's LANDSAT earth resources satellite is a new, powerful 
source of information that can contribute substantially to more 
informed decision-making. The LANDSAT is continually circling 
the earth, returning to the same point every 18 days. It is pro- 
ducing multispectral data that can provide a user with information 
at low cost for large areas. This information can be in the areas 
of geology, land cover, land use, water quality, crop conditions, 
etc. 


Other relatively new geographic computer analysis technologies 
have been emerging as a substantial aid in utilizing data from 
satellites and the more standard air photos, maps, and statistical 


^Myers, R, , 1976. "Data Needs in Missouri" - a presentation given at 
the conference on Future Directions for Earth Observation Data Manage- 
ment Systems, Washington University, St. Louis, MO, April 1976. 

^Council of State Governments, 1976. Natural Resources Data Needs 
Recommendations , Lexington, KY. 

Council of State Governments, 1975. An Evaluation of Natural Re- 
sources Data Products by State Data Users , Lexington, KY. 

Federation of Rocky Mountain States, 1973. Statewide Geographic 
Information Systems to Meet State Specific Needs, Denver, CO. 
Federation of Rocky Mountain States, 1972. Data Applications in 
State Land Use Planning, Denver, CO. 
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reports . ^ 

Many of these techniques are developed by the federal govern- 
ment and its private contractors and are not originally designed 
for operational use by state and local decision-makers. 

The demonstrated multisource data compositing approach of 
this project is designed for complex resource management and plan- 
ning problems such as shoivn in Figure 2.1.1. It is observed among 
states and counties in this region that multilevel decisions on 
complex problems may become bogged down in the "information prob- 
lem," and that the weaker the information, the more difficult the 
ultimate decision. This is particularly noticeable where several 
governmental jurisdictions are involved. The system approach de- 
signed by this project would improve interagency data- sharing and 
make possible mutual analyses such as agreement on the components 
of formulation of EIS procedures, systematic application of the data 
bank to state evaluation of local projects under various state 
laws, multiagency project evaluation under A-9S review procedures, 
and so on. 

Optimal use of CCT's calls for supplemental data.' The pro- 
ject found that while LANDSAT CCT input constituted the prime in- 
put to a regional data bank, it is really only con 5 )lementary to all 
other data and dependent upon other sources for complete answers to 
users. IVhile it can detect "reflectance characteristics" of many 
types land cover, it cannot account for "activity categories" which 
are essential to many analystic operations. In short, as the basic 
ground- cover description improves, so does the need for multi topic 
analysis in most kinds of functional planning. For these reasons, 
LANDSAT' s best utilization ivould seem to lie in raultisource regional 
data banks. 

In terras of "land use" identification, this leads to dual or 
multiple classifications containing both the "visible" (from remote 
sensing) and the "invisible," such as all categories urban zoning, 
economic activities, social data, legal boundaries, etc. 

Still another reason, for "data mixing" in a regional data 
bank is the heterogenous nature of original data coming from various 
agencies. The clientele of a regional data bank will include numer- 
ous bodies, providing peculiar forms of input and seeking numerous 
forms of output. 


“^The Federation has been involved in many efforts over the past 
seven years to improve states' capabilities to make decisions, 
using these computer mapping and satellite remote sensing technologies; 
for example: 

Information for Land Use Decisionmaking in Wyoming: A Computer 

Mapping Methodology (1974) . 

A Computer Mapping System for Regional Resources, Basic Industry 
and Public Facilities — Analysis and Projection. (1972). 

Jefferson County Open Space Study -- A Summary. (1975) . 
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FIGURE 2 . 1.1 
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industrial hydro geo graphic models of 
steam basins. 

Film product 
(1 250,000) 
Classified 
maps 

(1-24,000) 

Black ^ white 
medium alti- 
tude 

coverage 

Hydrology 
Ground water 

Various facility 
engineering maps 

Snow basin participation 
Forest conditions 

Stream flow networks 
Potential engineering 
Etc. 

Regional and urban studies, land 
use trends, urban corridors, in- 
dustrial locations, potential 
growth areas, etc. 

Film product 
(1:250,000) 
Classified 
maps 

(1:24,000) 

Black 5 white 
for detailed 
urban areas 

Land use surveys 
New subdivision 
5 zoning pro- 
posals. New 
urban facilities 
plans. 

Past land use patterns 
Present patterns (change) 
Imminent growth following water- 
sewer areas 

Political jurisdictions 
Land use plans, constraints 
Etc, 

Major recreation resources and 
development constraints, eacn 
season optimal human occupancy 
etc. 

Film product 
(1:250,000) 
Classified 
maps 

(1:48,000) 


Seasonal hximan 

occupancy 

sampling. 

Wildlife pattern maps 
Ecology map series 
Seasonal sports 
Land ownership 
Zoning constraints 

Critical land and water 
competition areas defined by 
states and localities (example 
Colorado Law 1041) for de- 
tailed review of new developments, 
mining subdivision, water 
diversions, etc. 

Film product 
(1.250,000) 
Classified 
maps 

(1:48,000 
for criti- 
cal areas) 

Black 5 white 
photography 

Subsurface min- 
eral/ground water 
maps. Industrial 
resources . 

Conservable resources (series) 
Potential industrial land 
and water uses 

Existing and potential urban 
or industrial uses 

Wide regional modeling of broad 
land uses and socio-economic 
distributions, future patterns 
expressing various planning 
scenarios, state development 
policies, federal coal,- shale, 
leasing and water development, 
etc. 

Film product 
(1.250,000) 


Population and 
employment, maps 
Industrial classi- 
fication maps 
Major resource and 
infrastructure 
existing 

Hypothetical devel- 
opment and popula- 
tion distribution 

Population density of various 
land uses 

Employment densities 
Potential employment growth, 
resources and infrastructure 
Environmental constraints 
Public facilities future pattern 
Zoning and density standards 
Etc. 

Up to date environmental 
impact statements, prepara- 
tion by private and public 
developers,, immediate 

Need to uti- 
lise above 
LAXDSAT 
products 

Need to uti- 
lise above 
airborne 
products 

Need to utilize 
above ground 
surveys 

Need to utilize above composite 
maps 


data surveys and compxlatLon 
of available secondary data 
showing plans and models of 
above activities, etc. 
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In these circumstances, the project. arrived at the 
cellular map, based upon the LANDSAT picture elements of 1,1 
acres. This defined the fine grid to which most other forms of 
data could be adapted, still permitting the nesting into larger 
cell sizes if needed for handling coarse regional data (economic, 
social, climate, etc.). There were other advantages in adopting 
the grid system. 

2.1.2 Recommendations 

Institutional Mechanisms . A principal question is, then, 
how to further enhance state and local governments' ability to 
utilize these new technologies on a day-to-day basis. This ques- 
tion of public technology, stimulating implementation of techniques 
(especially new, federally produced) by state and local governments, 
in many cases is a problem of developing adequate intergovernmental 
relationships and institutions resulting in transfer of a given 
technology. Roessner suggests that: 

"Public technology policy issues, then, rest at the 
nexus of two broad streams of public concern: the 

perceived need to stimulate the use of research and 
new technology by regional, state, and local govern- 
ments, and the need to assist these same governments 
as they assume greater responsibilities for the solu- 
tion of increasingly complex public problems. The 
problem for federal policymakers is how to develop 
and implement programs that will at once increase 
the availability and use of new technologies by 
state and local government and ensure that doing so 
actually helps solve the complex problems they face ."^ 

In the years ahead, the primary concerns of state govern- 
ment, in both the executive and legislative branches, are going to 
center on policy management — that is, enhancing their capacity 
for goal setting, evaluation and assessment, interinstitutional 
relations, and inter) urisdictional relations. A tremendous amount 
of knowledge and experience and good social theory are available 
to guide that process. Nevertheless, ways and means have to be 
found to get the leadership and the front-end resources together to 
engage in the massive process of institutional design and development 
if the great investments in research and technology made by the 
federal government are going to be used effectively to meet national 
needs. 


Itfhat IS lacking is the capacity to use the existing research 
and technology -- which in almost every case involves (1) creating 
new institutional capabilities (and especially those that will en- 
hance the ability of state and local governments and their elected 
leadership to participate more fully in the federal R^D agenda- 
setting and utilization process); and, (2) strengthening policy 
management capacity, including science and technology manpower 

^Roessner, J.D., 1976. Federal Policy and the Transfer of Technology 
to State and Local Government. 
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capabilities at the state and local level. 

There is a need to invest money and talent in the social 
architecture of research and technology in the design and 
development of better institutional mechanisms to use this 
existing research and technology. Then these same mechanisms 
v^ill help establish and sustain the communication between the 
different elements of the research and technology communnity — 
particularly research providers and the R§T resource users -- 
and this vi?ill, in turn, likely result in better research and 
operational useage of the technology. 

Data Banking . This project addressed the needs of "opera- 
tional" data banking and analytic services for "regional" programs 
of federal agencies, states, and multi county bodies. 

• An essential requirement is a broad range of data 
sources — physical, economic, and social — for 
regional development studies, wildlife studies, 
water quality and quantity planning, identifying 
critical areas of land conflict, and multifactor 
studies in transportation, energy development, etc. 

• Another requirement is ability to handle a broad range 
of formats — LANDSAT digital data; airborne sensing; 
freeform maps such as soils, crops, timber, or wild- 
life; human population characteristics; land oivnership 
patterns; local climatic and similar data types. 

• Ingestion of any input form or scale and converti- 
bility to any other scale. 

• Category transformation into different sets of symbols, 
different ordinal ranking, or various summarized cate- 
gories, etc. 

• Compositing ability to overlay numerous topic maps 
under given relationships, such as arithmetic, with 
different topic weights, or logical relationships, 
etc. 

• Or conversely, ability to analyze map relationships 
statistically, by subjecting digitized (preferably 

in cellular array format) maps to analysis, or multiple 
regression, etc., as elected by the analyst. 

While the original plan was to meet average resolution re- 
quirements with a grid of 2.5 acre cells, aggregative to 10 or 40 
acres, the LANDSAT processing subcontractor, Colorado State Univ- 
ersity, offered a fine grid of individual picture elements or 1.1 
acre cells. 
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Advantages : This had the advantage of matching the 

standard 1:24,000 USGS base map fin the LANDSAT cellular 
printout of 8x10 cells per inch), a convenient grid for 
other non- LANDSAT data compositing. Given the aggre- 
gative possibilities of the 1.1 acre cells to 10 acre 
cells (3x3x1. 1), this would be ideal for a regional data 
bank. 

Disadvantages : There were drawbacks in imposing such fine 

detail on LANDSAT training site work and verification. It 
was not necessary to use such a fine cell pattern for most 
subjects in regional planning. A significant factor is 
the LUDA program of USGS, i^rhich will be supplying 10-acre 
and 40-acre polygonal resolution elements (for Level II 
and Level I categories, respectively) as a major input to 
any regional data bank. 

In retrospect, any state centralized data bank should plan 
an "optimum" grid scale for extensive regional mapping, reserving a 
finer, nested grid for intensive local mapping — probably 40 acres 
and 1.1 acres respectively. 

To achieve truly operational status, cover large regions, 
and represent various agencies, a system would have to balance 
various parameters such as grid scale(s), land-use classification(s) , 
and the multi date survey calendar. It would have to do this with 
a view to the most common resource management problems at hand. 

Such balancing could not be attempted in this relatively 
brief project. As matters turned out, the LANDSAT Computer Com- 
patible Tapes picture elements were ambitiously adopted as 
the basis for a 1.1 acre grid system and also used for encoding 
all other types of data and experimenting with composite maps. At 
the conclusions of the project, some states reported that this cell 
resolution was too detailed and resulted in location registration 
problems for training sites and verification. This, together with 
the future available LUDA polygonal resolution of 10 to 40 acres, 
leaves a clear implication for aggregating LANDSAT pixels in the 
interpretation process and working toward LANDSAT output maps having 
10 ,or more acres per cell. 

Another problem in balancing the parameters concerns the de- 
tailed breakdown of categories. A great variety of land-use classi- 
fication was tried by the state teams, ranging from the original 
19 categories selected by the interstate workshop of April 19, 1975, 
up to 81 categories for one of the states. The various state lists 
were not compatible except at gross levels. This experience sho\?ed 
that it would be necessary for two or three adjoining states to begin a 
common Level I and Level II system from the outset, while still 
pursuing different categories at Level III. 
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Regional Assistance . An increased ability of the public 
sector — principally state and local governments — to use satellite 
remote sensing and related conputer mapping technologies on a con- 
tinuing basis to make decisions is required. 

This goal could be accomplished through aggregation of the 
user community (market) and transfer of the technology to this 
market by: 

• improving user (existing and potential) awareness of 
the technologies and their applications: 

• building user capabilities to utilize the technologies; 

• mobilizing existing resources in the region to provide 
maximum opportunities for technology transfer. 

Central to achieving these objectives set forth is the estab- 
lishment of a regional Earth Resources Technology Applications 
Assistance Center (ERTAAC) and subsequent development of operational 
programs. Such regional centers were recommended by a recent National 
Conference of State Legislatures report:^ 

A center would be widely used by various groups active in 
the Nest (see Figure 2.1.2). These groups would, in many cases, 
also provide resources for a networking effect. The center vvould 
be the hub — the catalyst — of a regional network of users, re- 
searchers, and suppliers. 

The focus of a Center's activities initially should be the 
applications of proven LANDSAT capabilities and related multisource 
computer mapping technologies. 


Functions of the center should include: 

1. User awareness and education — 

• formal and informal training (technical and 
applications); 

• information about the availability, costs, etc. 
of data and services; 

• availability of education resources such as 
speakers, audio/ visuals, etc.J 

• a central library for technicians and users. 

2. Information exchange and netivorking — 

• regular newsletter about innovations, applications, 
projects, activities, etc.; 

• a publications system providing reports, summaries, 
bibliographies, etc.;. 

• data exchange services with other information 
centers/ systems ; 


^ NCSL 1976, op. cit. 
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• an annual meeting for users and technicians pro- 
viding a forum for information exchange and con- 
tact (through exhibits^ with commercial and other 
suppliers; and 

• information referral services with libraries, 
organizations in the region, etc. 

3. Technical assistance and capacity building — 

o conducting user needs evaluations, 

e performing services and demonstratons, 

• personnel exchanges, 

• providing user access to data and hardware, « 

» assistance in establishing user in-house capabilities, 

• standardization of products and procedures to enhance 
operational uses, and 

• establishment of standard ground truth test sites in 
individual states. 

4. Mobilization of resources within the region -- 

• keep center staff at a minimum and use existing 
expertise, 

9 provide assistance in travel costs, etc., 

9 establish training/ assistance teams throughout 
the region, 

« coordinate regional activities as much as possible, 

9 provide referral and information exchange services. 


2.2 Toward Centralized State Services 


2.2.1 Findings 


IVhile this project involved a mix of functional state agencies 
and "technical" organizations, most state reports indicate the need 
and desirability for increased efforts of centralizing state data 
management. Although the project did not attempt a general review 
of the internal organization in member states, it found that informa- 
tion service organization was in flux and tending to centralize. 

A recent study by the Federation of Rocky Mountain States examined 
the issues and the evolution of land use planning operations in the 
member states.^ 

Following were the several state agencies in each state which 
appeared to have the closest concern for the kind of data system here 
being demonstrated: 


7 

Lament, W. and M. B. Jones. Land Use Planning in the Rocky Mountain 
Region: Background and Status . Federation of Rocky Mountain States 
Regional Background Paper, Denver, CO 1976. 
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FIGURE 2.1.2 


EXAMPLE REGIONAL ASSISTAI^CE USERS 
AND PARTICIPANTS IN THE mST 

FEDERAL GOVERNMENT LOCAL GOVERNMENT 


Department of Commerce 

• Environmental Research Labs 

• Weather Service 


• Counties 

• Regions 

e Service Districts 

• Cities 


Department of Housing & Urban Development PUBLIC SECTOR 


Department of Interior 



Interest Groups 

e Bureau of Reclamation 


Citizen Organizations 

« Geological Survey 

• 

Public-at“Large 

• Bureau of Land Management 

e Park Service 

PRIVATE 

SECTOR 

• Fish St Wildlife Service 


e 

Utilities 

Department of Agriculture 


Planners ^ Engineers 



Resource Firms 

e Forest Service 

« 

Research 

. Soil Conservation Service 


REGIONAL 

GROUPS 

.Department of Defense 


0 

Policy Development 

• Corps of Engineers 

• 

Coordination 



Research 

Department of Transportation 

0 

Planning 

Department of Health, Education 

and Welfare 

UNIVERSITIES 


• Public Health Service r State 

« Private 


STATE GOVERMIENT 


• Governors 

• Legislatures 

• Planning Agencies 

• Resources Agencies 

• Development Agencies 

• Human Resources Agencies 
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ARIZONA 

COLORADO 

MONTANA 

NEW MEXICO 

UTAH 

IVYOMING 


Arizona Resource Information System in 
the Department of Revenue 

Division of Planning, Office of Economic 
Planning § Development 

Colorado Division of Planning, Department 
of Local Affairs 

Colorado Land Use Commission 

Division of Planning, Montana Department 
of Community Affairs 

Department of Natural Resources 

State Planning Office 

Department of Minerals and Natural Resources 

Office of City Planning Coordinator 

Department of Community Affairs 

Wyoming Land Use Commission 

Planning Division, Wyoming Department of 
Economic Planning & Development 


A multisource data bank of value to various state and local 
agencies requires central technical functions. In this project 
the organizational situation differed widely among states, 
but generally it was found that either natural resource or planning 
agencies [whichever became the Lead Agency) began to assume more 
centralized responsibilities for the technology required. The project 
contributed to this evolution. Sometimes state data banking services 
had already been legislatively mandated to serve local jurisdictions 
and turned out to be useful also for other state and federal agencies. 
At the same time, the accessibility of federal and state aerial 
photo coverage and LANDSAT services tended to bring about this state 
central i z ati on , 


Many states in the West are now actively moving toward opera- 
tional usage of computer mapping and LANDSAT data. There is, how- 
ever, a definite evolutionary process, ultimately arriving at the 
broadest, most flexible usage of these technologies in making policy 
and program decisions; no western state has yet reached this climax. 

More particularly, state situations could be described as 
evolving up a "ladder" of data operations [see Figure 2.2.1): 

• Stages I and II .- where there exists only "ad hoc" 
assembly of data by each agency for its own purpose, 
generally resulting in duplicative and expensive 
efforts by various operational agencies, and eventually 
executive or legislative directives toward centraliza- 
tion. 
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• Stage III - where a general purpose "map library" 
function develops in an agency having widespread 
responsibility for state planning or for some exten- 
sive repetitive survey of natural resources, etc., 
procurement of remote sensing, possibly also carto- 
graphic standardization, gathering and indexing of 
various agency data products, etc. 

• Stage IV ” iifhere analytic and modeling services are 
added to the foregoing map library, making possible 
more complex map analyses, regional change studies, 
and provision of information services to other agen- 
cies, to the governors and legislative committees, etc. 


FIGURE 2.2.1 







STEPS TO COMPREHENSIVE STATE DATA 
APPLICATIONS SYSTEMS 





multisource, multiagency 
use of information system 



muitisourcc sxngle agency 
system, serving single 

with analytic ^ervxces 


few sources single agency, 
ad hoc purposes 

purpose 


ad hoc data search 
by :> ingle agencies 




2.2.2 Recommendations 


' Information, dissemination, demonstrations, and technical 
assistance to states and localities — either unilaterally or via 
a regional approach — are required to move them up this evolutionary 
ladder, to build their technical capacity and to increase their in- 
terest in utilizing the technology. 

Continued upgrading of LANDSAT products and use in combina- 
tion with multisource data will encourage state efforts at standard- 
ization, coordination, and centralization of data management. 

Additional training of state personnel, applications demonstra- 
tions, standard procedures, improved accuracy and lowering of costs 
will also assist in developing state efforts. 
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2.3 Toward Operational Use of LANDSAT 


2.3.1 Findings 


As might be expected for a project of this type, the results 
were a mixture of good and bad. The most disappointing aspect of 
the project stems from our inability to achieve the maximum capa- 
bility from the LANDSAT imagery. This inability to maximize the 
results was caused by lack .of communications, inejq)erien'ce on the 
part of all project participants, by the lack of time and funds to 
correct deficiencies which were noted during the first processing 
of the data. Even with the most e^qserienced personnel, an iterative 
mode of operation is required to achieve the best results. 

On the positive side, the discrimination capability of multi- 
date LANDSAT imagery was very encouraging. Separation of forest 
species, residential types, vegetation density/ condition classes, 
and many other Level III types appears to be possible. 

Although there are technical questions which remain to be 
answered, the largest question is one of hoi-j to achieve operational 
status. LANDSAT and other remote sensing systems can provide useful 
information for land use planners and natural resource managers. 

We have yet to accomplish a mode of operation which is timely, cost 
effective and consistently reliable. 

All indications from this project and other projects of a 
similar nature show that LANDSAT has the capability of mapping land 
use wherever such use affects the land cover characteristics in a 
consistent manner. IVhen the land use does not affect land cover 
characteristics, LANDSAT is, of course, of no value. From the com- 
ments received from the state participants associated with this 
project, this appears to be a serious limitation from their point 
of view. In many instances, state and local organizations are more 
concerned with oi>mership and activity characteristics of land use 
than with those characteristics which affect land cover. In other 
words, they would like to know, -not that an area is covered with 
grass and trees, but whether it.'.is used as a golf course, cemetery, 
park, pasture, or some other activity. 

Clearly, LANDSAT promoters are faced with a need to communi- 
cate more effectively the potential benefits of LANDSAT data in 
order to obtain the support from state and local agencies which will 
be needed before an operational status will be achieved. This will 
be considered again under recommendations. 

All in all, the results from this project show that many land 
cover characteristics can be classified using LANDSAT data. The use 
of multidate classification enhances this capability under many cir- 
cumstances. Although we need to learn more about the temporal 
characteristics of vegetation type and land use (including activity 
types of land use) to maximize the potential from multidate process- 
ing, this project has shown that the capability is there. 
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It is important for all participants in this project to 
realize that the maximvun capability from multidate LANDSAT imagery 
was not achieved. In the one instance when the opportunity arose 
to retrace s.teps, obtain a better s"election of classes and better 
training data to represent them, very satisfactory results were 
achieved. This was accomplished for the Buffalo, Wyoming, quad- 
rangle after the preparation of the Wyoming State final report. 

In this one instance, there also was the opportunity to compare 
the multidate- results with single-date classification, using the 
best single date possible for classifying the land uses in the 
Buffalo quadrangle. The multidate classification was clearly 
superior. A quantitative evaluation of the superiority in terms 
of accuracy of classification has not been made, but it is apparent 
from a side-by-side comparison that the multidate map is more con- 
sistent with known land uses in the quadrangle. The single-date 
classification results themselves, however, appear to be quite 
satisfactory . 

2.3.2 Recommendations 


As previously stated, many of the efforts toward operational 
use of LANDSAT lie in developing nextf institutional mechanisms for 
multisource data management, oriented to "real" user problems. 

Products. Much work remains to be done before routine use 
of LANDSAT data will be accomplished in a timely and cost effective 
manner, with output products being enthusiastically received by the 
land use planner and natural resource manager. 

Perhaps as important as any change which could be made in 
the use of LANDSAT would be an in 5 >rovement in the output products 
to make them more palatable to the users. This could be accomplished 
in many ways, but perhaps one of the most important would be to pro- 
vide the user with an output product generated not only from LAND- 
SAT information. If those organizations generating LANDSAT products 
would combine them with digitized information relative to land owner- 
ship, zoning, and activity types of information, the final product 
provided to the users might be much more acceptable. 

Training Site Procedures . It is recommended that a training 
site manual be written in a standardized form for use by nontechnical 
personnel, concerning the following subjects: 

• Selection of land use and cover categories in considera- 
tion of the end purposes of a study and availability of 
other more accurate information on functions and activi- 
ties on the land to be used in a compositing system. 

• Particularly, selection of LANDSAT categories for their 
spectral characteristics. 
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• Hierarchiacal classifications of land uses, cover, 
and functions, permitting appropriate aggregation 
in the LMS and in the composite mapping programs. 

• Rules concerning homogeneity of land cover conditions, 
mixture of vegetation, etc., based on lessons learned 
in this and other lANDSAT exercises. 

‘ j- 

• New rules concerning multi date ground/ vegetation con- 
ditions . ;i 

« Rules concerning the location and size of training fields, 
concerning amount of data needed, size of sites, number 
of sites, site-to-site variability for a given class, 
slope/aspect considerations, etc. 

© Rules on dociimentation of training data: field mapping 

and documenting procedures, particularly seasonal 
tracing of changes; standard annotation form for 
description of all pertinent field characteristics; 
and optimum use of other sources of data, aerial photog- 
raphy, etc. 

Standard Training Sites . In order to approach a fully opera- 
ting system, for optimum use of LANDSAT and an all-source data bank 
over wide regions, the root problem of ground truth must be solved. 
Many weaknesses were experienced in organization of this work on a 
first- time basis ^^ith inexperienced personnel, which will reoccur 
ivith each new LANDSAT application. Therefore, a key recommendation 
for follow-on work is to establish a "standard test site system" vi?ith 
numerous permanent ground truth survey sites under repeated observa- 
tion for changing land uses. This would go far toward practical 
LANDSAT use in this large region of public and private lands. A 
system of standard training sites will provide first grade seasonal 
or annual coverage of the changeable land uses, cover types, water 
and moisture status, etc. Sites ifould be distributed in strategic 
locations within each LANDSAT scene. Any new user of LANDSAT need 
only call upon the new bank of standard training -site data for the 
ground truth for practically any LANDSAT or remote sensing survey. 

More particularly, the reasons and advantages of such a 
program are: 

9 The optimal ground resolution of LANDSAT is relatively 
coarse but very cost effective for wider regional appli- 
cations . 

© Similarly, it is made increasingly effective by repeti- 
tive coverage. 
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• LANDSAT's relative freedom from geometric distortion 
is a large advantage in repetitive regional coverage. 

• LANDSAT's specia-1 ability to discern vegetation con- 
ditions, biomass, water conditions, etc. is ideal for 
seasonal and/or annual repeat coverage. 

• The cellular output format of LANDSAT CCT's makes 
possible a basic grid suitable also for supplemental 
mapping of various land use "functions," slope and 
aspect, etc. to optimize land use cover classification 
for a data bank. 

• The new LANDSAT systems, C and Follow-on, provide more 
and different reflectance bands and will require a great 
deal of new ground truth to match them. 

• Any inadequacies in LANDSAT discrimination may be augmented 
by compositing with other information on crop conditions 
and production levels, forest types and conditions, urban 
land uses, mining activity zones, construction zones, and 
recognizable geodetic reference points. 

Concerning the locations of standard sites, experience in this 
project suggests that standard sites be located in each LANDSAT scene — 
enough sites to represent" each important category. For example, if a 
state needed 16 LANDSAT frames for complete coverage, and if there 
were 20 land use and cover categories of interest, the required number 
of sites would be something less than 16x20, because some larger sites 
could be selected which have subareas that cleanly represent the cate- 
gories . 


In repetitive system, there would be certain economies. For 
exan^le, the first round of ground truth would have to survey every 
category of interest. It would establish (1) the permanent character- 
istics of slopes, aspects, soils, major forests, riverine vegetation, 
etc., and the changeable characteristics of agriculture, water 
status, surface disturbance from mining, construction, etc. However, 
future repetitive surveys could focus only on the changeable land uses 
and vegetation, at a large savings in work. 

Such a system would permit the use of best technical methods — 
usually also the most economic, if used repetitively. Various agencies 
would become adapted to the system, contribute their survey methods and 
data, and call upon the central system to support any LANDSAT or other 
remote sensing project they wish to initiate. 
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Essential technical dimensions of the system might include: 

• A sufficient number of designated sites in each LANDSAT 
scene. 

• Areas of 160 acres to several square miles each. 

• Grouping of the test sites for economy of field work 
ivhile preserving the homogeneity of individual cate- 
gories of subareas. 

• Most characteristics of the categories located in the 
land use/cover areas . 

• Location in characteristic land forms, slopes, aspects, 
soil zones. 

• Possible categories — 

3 - 6 rangeland types 

5-10 crop types 
3 - 6 forest types 

3 - 6 local water "states" 

2 - 3 types of land surface disturbance 

slope, aspect, general elevation 
soil classes, moisture conditions 
surficial geology 

crop and forest condition indicators, disease, biomass, 
leaf moisture, erosion state, etc. 

Starting such a system would require the appointment in each 
state of a "permanent survey agency" v^ith a standard procedure manual, 
common training, close working relationships with the remote sensing 
processing facilities. 

The system would also require one or several remote sensing 
processing facilities, under appropriate interstate and state-federal 
arrangements. In the beginning, a single center of strategic planning 
and training would be needed to ensure uniformity of approach. 


2.4 The LANDSAT Mapping System 


The LANDSAT Mapping System (LMS) was developed for this project 
by Colorado State University.^ The major and minor operations built in- 
to LMS are illustrated in the CSU portion of this report and will not 
be repeated here. Among its objectives were: (1) to simplify human 

interactions, (2) make possible efficient multidate LANDSAT applica- 
tions, (3) handle a variety of locations of training data, (4j efficiently 
develop class signatures, (5) perform selective masking of training fields 
and (6) focus statistical analysis on the remaining training data. 


^The LMS System was evolved from Colorado State University's earlier 
Pattern Recognition Program (RECOG) by Dr. Lee Miller and Dr. Eugene 
Maxwell. 
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2.4.1 Findings 

Data Handling . The data-handling capabilities of LMS are 
not entirely satisfactory. It is time consuming to obtain given 
portions of data in the desired format for the different steps of 
LMS and still leaves too many opportunities for human errors. 

Even though this system has achieved an order of magnitude improve- 
ment over the old RECOG or LARSYS type system, it still falls short 
of the capability of interactive systems. CSU now has an inter- 
active computer capability and as time and funds permit, an improve- 
ment in the data handling capabilities will be sought. 

Geometric Corrections . The geometric corrections achieved 
by LMS appeared to be adequate for quadrangle size areas. If it 
were desired to correct and display an entire scene as a single 
image, it is doubtful that these corrections would be adequate. 

In general, however, the geometric correction capability of LMS 
is quite flexible and effective. 

Radiometric Corrections . The only radiometric correction 
which LMS had during this project was in the form of low pass 
filtering. This achieved some correction for radiometric banding, 
but it was not adequate to correct serious problems of this type. 
Future modifications of LMS should provide for banding correction 
and for removal of atmospheric effects. A recent modification of 
LMS, not available during the lifetime of this project, provides 
for the replacement of missing pixels or entire missing scan lines. 

Signature Development . LMS has excellent flexibility for 
computing training- data statistics and in modifying and correcting 
that data to achieve the best signature. The speed and efficiency 
of carrying out these operations, however, leave something to be 
desired. Again, a more truly interactive mode of operation with 
LMS might correct this deficiency. 

Classification Algorithm. The classification algorithm 
of LMS is as effective and accurate as any other software system 
currently being used. It employs a maximum likelihood algorithm 
which could be easily modified to use a priori probabilities and 
other weighting factors to achieve the best accuracy possible. 

Being a software system, it is not as fast and cost effective as 
hardware systems, but this is probably its only serious deficiency. 

Display . The display capability of LMS, being limited to 
microfilm or line printer types of output, is deficient primarily 
in terms of product appeal. Technically, the LMS products are 
as accurate and useful as products supplied by any system, but we 
cannot ignore the appeal of photographic color displays, which LMS 
cannot produce. Incorporation of some sort of a color display capa- 
bility with the system is highly desirable. 

Costs, Because of the developmental aspects of using LMS 
during this project and the inexperienced project personnel, the 
costs discussed do not reflect an operational mode. Costs and 
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comparisons of LANDSAT versus use of other data types are 
discussed in each of the reports in the appendices. CSU 
reported computer processing costs for classification at 
$300 per quadrangle, or about $5 to $6 per square mile. 
Complete analysis costs, including man- time, would be much 
higher, but total costs are unknown' (and presumably lower) 
for an operational mode. 

2.4.2 Recommendations 


Class Selection . To achieve the most effective results 
from LANDSAT data, it is important that careful consideration be 
given to the requirements imposed by the assumption of a normal 
distribution of radiance values. This assiunption requires that 
the radiance detected by the satellite sensors should come from 
the same mixture of materials on the ground for each pixel with- 
in a given class. This, and considerations of the reflectance 
characteristics of materials, should be used as a guide in the 
selection of classes at the outset of any project. Problems 
associated with improper selection of classes are undoubtedly a 
major factor in the current popularity of unsupervised, cluster- 
ing modes of operation. The unsupervised, clustering mode of 
operation automatically takes into account the factors noted above. 

Total dependence on unsupervised clustering, however, will 
increase the problem of convincing users that LANDSAT can provide 
useful information for their applications. If there is difficulty 
in convincing state and local agencies that land cover of known 
categories are of value, vtfhat chance will there be to convince 
them that some unknown and changing set of land cover characteristics 
can be fit into their scheme of operations? Furthermore, the dis- 
tinct advantages of multidate classification could well be lost by 
an unsupervised, clustering type of analysis. 

The preparation of a class selection and training data manual 
might help to eliminate some of the problems ivhich stem from improper 
class selection. 

Training Data Selection . Similar to the problems of class 
selection, the selection of training data by inexperienced and un- 
trained personnel can detract greatly from the results achievable 
with LANDSAT data. The preparation of a manual, clearly defining 
the requirements for training data (to optimize final mapping results) 
IS highly desirable - 

Signature Development . Considerably more research needs to 
be done relative to signature development and the optimization of 
both the supervised and unsupervised algorithms. Although inter- 
active modes of operation, which depend on testing the classifica- 
tion results, are generally quite effective, we cannot depend on 
such methods if we are to achieve repeat results for essentially 
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the same classes at different times in different locations. Much 
research is under way at the present time relative to signature 
extension and signature problems in general; so there i-s reason to 
bel-ieve that many of these problems iirill be solved in the near 
future . 


Software Documentation . Documentation of the LMS software 
is an absolute requirement before it can be used effectively by 
other organizations. Adequate users' and operators' manuals are 
required to train new personnel as they begin use of the system. 

CSU will make this system available to anyone desiring to incorpor- 
ate it into their computer operations. This cannot be effectively 
accomplished, however, without adequate documentation. 



3.0 STUDY AREAS 


3.1 General Study Areas 

Figure 3.1 outlines the large study areas selected to represent 
characteristic natural resources and economic activity zones in the 
Rocky Mountain region. 

The Montana- Wyoming area contains some of the vast Powder River 
coal deposits. Energy developments will significantly affect the pre- 
vailing grazing and farming. This area, part of the Great Plains phys- 
iographic region, is mostly rolling grass and brush prairies with scattered, 
wooded uplands. 

The Utah area was selected because it is a rapidly urbanizing cor- 
ridor along the Wasatch mountain range. A major issue is the transi- 
tion of agricultural and grazing lands to urban and industrial uses. 

There is commercial timber in the Wasatch range and large wetlands a- 
long the Great Salt Lake. The capitol city. Salt Lake City, is in the 
study area. 

The Colorado-New Mexico area is characteristic of southern Rocky 
Mountain physiography. It covers part of the Rio Grande river valley. 
Grazing, irrigated agriculture and forestry are important land uses. 

Second home developments are increasing in the study area. Santa Fe, 

New Mexico's capitol city, and a number of other small population centers 
lie in the study area. 

The Arizona study area covers the Phoenix metropolitan area, where 
widely spreading urbanization is moving into areas of irrigated agri- 
culture and desert. Located in the Basin and Range physiographic region, 
the area has a hot desert climate. 
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Within these large study areas, each state selected four smaller 
test sites. These served as foci of the detailed data collection and 
analysis. 

3.2 Specific Test Sites 

Within the large characteristic study areas, each state selected 
four U.S. Geological Survey 7^2-minute quadrangles (scale 1:24,000) for 
detailed training site and computer analysis work (Figures 3.2 through 
3.7). 


Selection of these test sites was accomplished early in the project 
so that best cloud-free multidate LANDSAT and high altitude aerial data 
could be sought. Numerous training sites were then selected in these 
quadrangles. Additional information about these sites may be found in 
the individual state reports. Appendices C through H. 

Each state selected one of their test sites (shaded in the following 
figures) for conducting the multisource computer compositing described in 
Section 6.0. 

3.2.1 Arizona Sites 


The Arizona test sites are all located within the Phoenix 
metropolitan area. This area is undergoing rapid urbanization of 
agricultural and open-space lands. The test sites offer a mix of 
land forms: steep, rocky mountains, rolling hills, flats, and 

watercourses. The hot desert climate and resultant ecosystems 
offer unique constraints for agriculture and urban development. 

Citrus crops are the most significant agricultural activity. 

The Tempe site is highly urbanized with scattered agricultural 
lands. Most of the site is flat. There is a group of steep mountains 
in the center of the site, and the Salt River -- the major watercourse 
in the region -- bisects the area. 

Paradise Valley is also an urban area and includes part of 
the cities of Phoenix and Scottsdale. The several scattered mountains 
are surrounded by residential land uses. 

The Tolleson site is less urbanized, being mostly scattered 
small communities and extensive agriculture. It consists of almost 
completely flat, mountainless terrain, with one major watercourse 
in the southern-most sections of the site. 

The comprehensive test site, Hedgpeth Hills , is mostly rural 
in nature with subdivisions encroaching on agricultural lands. 

Several streams flow through the area and steep mountains and 
ridges persist in the northern sections. 
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3.2.2 Colorado Test Sites 


All sites are situated in the San Luis Valley of Southern 
Colorado. The sites are representative of land uses and vegeta- 
tion reflecting the entire study area. 

Zapata Ranch characterizes the tightly stacked vegetation 
zones on the west slope rise of the Sangre de Cristo Mountains. It 
also provides instances of rangeland and mass wasting processes as- 
sociated with the sand accumulations in that sector of the valley. 

Alamosa West contains commercial and semi-urban land uses, 
some agriculture, and a sample of the greasewood-saltmarsh complex 
which is scattered throughout the valley on the less well-drained 
portions. 

Manas sa quadrangle exhibits examples of riparian and basalt 
butte communities, and includes the best cross-section agricultural 
and small-town land uses. 

The Fox Creek quadrangle, selected as the comprehensive test 
site, is situated on the western edge of the San Luis Valley. It 
is located partly in the San Juan Mountains. Here the mountains 
rise in a very moderate ups lope from east to west, capped by a re- 
sistant basalt layer. Cutting through this flat plateau top is 
the Conejos River and several minor drainages. As is the case with 
the Zapata Ranch quadrangle, agricultural lands are few, scattered 
along the Conejos River. The predominant vegetation type, adapted 
to the lower elevations, is open, dry range characterized by bunch 
grasses and small shrubs. Mixed with this range type are stands of 
pinyon- juniper, Ponderosa pine, and some sagebrush. In the higher 
valleys, Douglas fir/White fir communities occupy steep northern 
exposures. Aspen groves and meadow lands are found on higher ele- 
vation sites in the drainages. The high tabletop land of this up- 
slope is dominated by mixtures of spruce, fir, and aspen. These 
areas were heavily logged in the early 1900 's which accounts for a 
variety of vegetation patterns, and it is the problem of revegeta- 
tion which provides the compositing exercise of revegetation index 
mapping. 


3.2.3 Montana Test Sites 


The test sites in Montana are underlaid by Powder River coal. 
Surface mining and large-scale energy developments are underway with 
more planned for the future. 

The Beaver Creek School site's major land uses are dryland 
farming on the valley floors. Timber is produced in the shallow 
soils of the deeply incised drainages and is an important economic 
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crop to the area, as is hay production on low terraces and 
floodplains. Grazing land dominates the moderately steep 
rolling uplands. 

On Poker Jim Butte , timber production and grazing of 
native rangeland are the principal land uses. Small amounts 
of hay are being grown in the floodplains. The timber is con- 
fined to the steeper slopes but is quite productive and impor- 
tant economically. Recreational use will drastically increase 
with energy development. 

On the Decker site, rangeland dominates the quadrangle with 
very limited amounts of dry cropland. Some irrigated hay is grown 
along the lower fens, terraces and floodplains of the Tongue River. 
Intensive energy development is currently underway with the Decker 
mine expanding its production and applying for new mining permits 
on the east side of the Tongue River. The Tongue River Reservoir is 
important for water resource management, irrigation, recreation, 
and potential energy conversion uses. 

The major land uses of Colstrip Southeast , the comprehensive 
test site, are rangeland with limited dry crop land and timber pro- 
duction. Intensive surface mining of coal is being undertaken 
and several levels of development of the power plant at Colstrip 
are under construction. A variety of stages of reclamation of the 
surface-mined areas is present. Interest in the quadrangle stems 
from energy development and conversion and its potential impacts 
on land use, recreational opportunities and the quality of the 
surrounding environment. 

3.2.4 New Mexico Test Sites 


Characteristically, the mountains in the test sites are 
forested and valleys are brush covered with a smattering of pinyon 
and juniper. A scattered population inhabits this rural area and only 
in the Santa Fe quadrangle, which contains the city of Santa Fe, is 
there a large concentration of people. Valley agriculture and graz- 
ing dominate the region with raining and forestry being important 
locally. Land ownership or management responsibility is under the 
control of the U.S. Forest Service, Bureau of Land Management, several 
northern Indian pueblos, and private individuals. In selecting the 
quadrangles, an effort was made to focus on one of several charac- 
teristics within the region. 
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land uses include agriculture and commercial forestry, 
Questa and Guadalupe Mountain are adjacent to one another and 
represent extensive sage grazing lands, a large molybdenum open- 
pit mine with related mill and settling ponds, as well as diverse 
commercial timber stands. The Bureau of Land Management adminis- 
ters much of the sage grazing land. 

The Santa Fe quadrangle, the comprehensive test site, is 
an urban area at the base of the Sangre de Cristo Mountains and 
contains the city of Santa Fe, Strong population and development 
pressures are being applied to an area of limited resources. The 
city lies within, and to the south of^ the Santa Fe River flood 
plain, a gently sloping valley. Topographic relief is great; with- 
in 10 miles the elevation drops from 8,500 feet in the east to 6,500 
feet in the southwest. 

Water is primarily supplied through wells but an increasing 
emphasis is being placed on reservoirs along the Santa Fe River. 

Many of the soils have high shrink and swell characteristics and 
low permeability which impairs the development of an area without 
help from city utilities and proper foundation development. Develop- 
ment over arroyos increases the flood hazard problem, and efforts to 
combat such development are being built into city zoning guidelines. 

3.2,5 Utah Test Sites 

The Utah test sites are in the Wasatch Front, an area of 
dynamic urbanization and dense population. 

The Tremonton quadrangle represents something of the anti- 
thesis of The Wasatch Front urbanizing trend. Farms there are 
larger, the agricultural economy is prevalent, and the distance 
from urban clusters is great enough to reduce the threat of rapid 
change. 

"The Dromedary Peak quadrangle is the only representative of 
mountain environments. It is typically mountainous northeastern 
Utah, with steep slopes, many aspects, and covered by a mixture of 
conifers, deciduous trees, brush, meadow, and bare rock. The Salt 
Lake City South site is an ideal example of rapid urbanization and 
small field agricultural activity as well as land-use change near 
the city core. 

Farmington , as the comprehensive test quadrangle, represents 
the epitome of dynamic variety. It includes an arm of Great Salt 
Lake, which itself is undergoing constant fluctuations in water level 
and economic/political battles for industrial, recreational, wildlife. 
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and hydrologic demands. The transition from irrigated agriculture 
to urban/industrial/transportation corridor uses is of great concern 
to planners in the area for three reasons: (1) loss of farmland; 

(2) scattered and uncontrolled development; and (3) encroachment into 
sensitive foothill areas threatened by flooding, faulting, and slump- 
ing. Foothill development further compounds the hazards of instabi- 

threatening land and structures and seriously reducing big 
game winter habitat. 

3.2.6 Wyoming Test Sites 

The Wyoming test sites offer a variety of land uses, including 
energy (coal) development, urban development, water storage, forestry 
and grazing. 

quadrangle is in gently rolling terrain, underlaid by 
thick coal deposits. It is the site of an active strip coal mine, 
one of the few in the Wyoming Project area. Riparian vegetation 
and agricultural land uses, as well as grazing areas, are present 
on the quadrangle. 

Hunter Mesa is close to the Buffalo quadrangle, but is in 
the lower reaches of the Bighorn Mountains. It is wholely within 
Bighorn National Forest, is mostly wooded (lodgepole pine, ponderosa 
pine, aspen), and contains some significant mountain meadows or 
"parks." Pressure for new recreational facilities on sites contained 
on that quadrangle is likely, as population increases in Buffalo to 
the east. 

Lake de Smet, West , is to the north and west of nearby Buffalo. 
This quadrangle covers much of Lake de Smet, a natural reservoir not 
far from the base of the Bighorns. It is in disected table land, with 
low to moderate relief. The most extensive land use is low-intensity 
grazing, although irrigated cultivation does occur along some of the 
creeks that cross the quadrangle. It is likely that coal mining will 
soon take place along the margins of Lake de Smet. 

Buffalo quadrangle is the comprehensive test site and is in 
gently rolling terrain at the base of the Bighorn Mountains. It is 
the site of Buffalo, Wyoming, one of two cities in the test area 
whose population exceeds 500 (Buffalo's 1970 population was 3,394). 
Strippable coal deposits underlie portions of the quadrangle and 
coal mining is presently occurring some six miles north of Buffalo. 

As regional development proceeds (primarily because of increase in 
coal extraction and related technologies), it is expected that 
Buffalo will be heavily impacted. It was primarily for this reason 
that the Buffalo quadrangle was selected for inclusion in the study. 
IVhile expansion of the town of Buffalo seems inevitable, there are 
many constraints to growth which lead to vexing planning problems. 

Land uses in the quadrangle include residential, commercial agri- 
culture, and grazing, as well as some riparian vegetation. 
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On key questions of policy and technique, the Federation Conmiittee on 
Earth Resources Technology Applications (representing scientific and administra- 
tive bodies in the region) periodically met and reviewed the project. 


4.2 Geographic Approach 

Two -representative interstate resource areas of approximately 5,000 square 
miles each were selected, particularly to cross state boundaries and link four 
of the states to encourage use of common remote sensing procedures, land/cover 
categories and analystic procedures. The Montana and Wyoming test region 
included active surface mining coal lands, irrigated and -dry farming and 
limited watercourses subject to intensive water engineering. The Colorado 
and New Mexico test region followed the main drainage of the Rio Grande 
River. The two intrastate test regions of Utah and Arizona were defined as 
dynamic metropolitan areas where conflicts would continue to occur between 
existing and potential land and water uses. 

Within each state's large test regions, four quadrangles were selected 
as having characteristic types of land use and vegetation. These quadrangles 
were used for training site work and for LANDSAT land use/cover mapping. 
Ultimately one of the four quadrangles became the target area (comprehensive 
test site) for full land use analysis. Correspondingly, final LANDSAT veri- 
fication work was conducted in these quadrangles (see also Chapter 3 of this 
report) , 


4.3 Procedural Approach 

Four general phases were followed as outlined in the following svimraary 
of the project procedurec 

1. Data procurement and preparation for later analyses. 

• Large study areas were reviewed and four 7%-minute quadrangles 
selected in each state. 


Initial land use and cover classes were selected, beginning with 
classes of Level I adapted from the USGS Anderson-Hardy table. ^ 
These were later modified to those shown in Figure 4.3. 


iy 


•- Initial rules and procedures were defined for selecting sites and 
describing and documenting site data, 

• High altitude color infrared photography was procured from NASA 
to aid in the training-site work. 

• Multidate LANDSAT imagery was selected for three or four seasons 
(1974) for each test area. 


0 Compositing exercises were formulated for the comprehensive test 
quadrangle in each state according to a characteristic resource 
management problem existing in the area. 


Anderson, J.R., et. al., A Land Use Classificat ion System for Use with Remote 
Sensor Data . U.S. Geological Survey Circular 671, 1972. ^ ^ ~~ ’ 
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© Corresponding types of multisource data (maps) were selected 
and converted into 1.1 acre grid cells, with various weights 
assigned for relative importance. LANDSAT cellular maps were 
edited for selection of certain features to be used in the 
composites. 


FIGURE 4.3 
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2. LANDSAT data analysis. 

• A new LANDSAT Mapping System (LMS) for CCT processing was 
developed by Colorado State University. 

® Initial grey scale LANDSAT printouts at 1.1 acre (pixel or 
picture element size) were produced for state quadrangles, 
showing unsupervised variations in reflectance. 

• Classified land use/cover maps were produced for several 
quadrangles in each state, using their given third and fourth 
level land- use categories. 

• Field verification work was conducted by the states, using 
common forms and procedures as adapted in regional workshops. 

3. Multi-source computer map compositing. 

• Computer work was performed by Los Alamos Scientific Labs 

for four of the states, using a special interactive grid mapping 
program. Generalized Map Analysis Planning System (GMAPS) , 
developed from the Composite Mapping System (CMS-II) by Dr. 

Keith Turner, Chairman of the FRMS Advisory Committee. 

• Two states — Utah and Montana — performed their own com- 
posite mapping with original programs (Resource Analysis 
Program (RAP), and Environmental Resource Geo- Information 
System (ERG IS) . 

4. Evaluation and reporting. 

• Three formal advisory meetings were held to discuss and 
evaluate project progress and findings. 

• Several other review sessions and state workshops were held. 

• Individual reports were developed by each of the six states, 

CSU, and LASL (see appendices) . 

• This final report will be given wide distribution to 
technicians and policymakers alike. 


4.4 The Multisource Approach 


One of the project objectives, of course, was to utilize LANDSAT data 
in conjunction with data from other sources to demonstrate procedures and 
applications. Several factors warrant discussion of this approach. 

The decision to use the finest grid size of 1.1 acre cells followed the 
initiative of Colorado State University in programming a new LANDSAT Mapping 
System for the purpose of efficient output of pixel^ level maps (see the CSU 
report. Appendix A) . This cell size was finer than necessary for most of the 


^LANDSAT picture element, also about 1.1 acres. 
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other relevant data to be used in compositing, yet permitted the preservation 
in the data bank of full details for other purposes. The 1.1 acre cellular 
format provides the most discrete resolution for LANDSAT mapping, and offers a 
convenient line printer output scale of 1:24,000 which corresponds with the 
7%-minute quadrangle map series of USGS and, therefore, enables use of over- 
lays. A very significant commitment was that of using multidate (seasonal) 
imagery to demonstrate the efficiency of the new LMS software and provide as 
complete a series of land cover signatures as possible. This began 'W-ith a 
search for good LANDSAT coverage during three or four seasons of on'e year, 
and the resulting adoption of an optimum sequence during 1974 (three seasons 
in some states and four in others) . 

Several workshops were held in which training-site analysis procedures 
were formulated. Subsequently, forms and instructions for documentation were 
developed by CSU and the Federation and published in the periodical project 
journal. The Remote Sensor . 

The following types of training-site problems were anticipated: 

1. Avoiding or fully describing variability within classes. 

2. Documenting any significant slope, aspect and ground moisture 
variations. 

3. Variable tree crom density and variable soil backgrounds. 

4. Minimum required site sizes and number of sites. 

5. Location and documentation of slope and aspect. 

6. Possible changes of land use since the 1974 LANDSAT image 
dates and how to research the original state. 

7. Locational referencing problems. 

The compositing objectives required the selection of other non- LANDSAT 
topic maps — social, economic and physical. These formulations were de- 
veloped through several project workshops and additional meetings of interested 
agencies in each state, under the guidance of the state lead agencies. Com- 
positing exercises were developed to analyze relevant resource planning problems 
in each test site area. (This is further discussed in Section 6.0.) 

In this procedure, the land-us e/ cover information from LANDSAT served 
variously as positive or negative overlays or as land-use constraints in 
relation to the geological or urban characteristics, etc. shown by the other 
maps. Studies have been made in each state and in each federal agency con- 
cerning land-use classification systems for various purposes. Every state has 
several agencies concerned with this. Consolidation studies have been made by 
the Regional Planning Council of the Federation of Rocky Mountain States. _ But 
the user diversity persists, and the more appropriate agreement would be in a 
process for interchanging areh land-use studies -and observing essential nesting 


%RMS Land Use Classification Project. 1971. 

FRMS Toward a Common System of First Order Land Use Classification . 1971 
FRMS Toward a Common System of Land Use Classification . 1972. 

FRMS Data Applications in State Land Use Planning. 1972. 
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definitions (levels of classification) plus establishing common geographic 
boundaries for polygon- grid and grid-to-grid nesting. In each functional 
area — agriculture, natural resources, etc. — there will ineva-t-ably occu3 
different classi-f-ications of particular activities or conditions of re- 
sources. Therefore, a "general purpose" system must be approached from 
"elementary indicators" which can be mixed into more complex activities 
and functions . 

Remote sensing can generally see only elementary indicators of certain 
kinds. It is equally important to use low-altitude photography, geologica] 
data, water data, industrial and urban data and other kinds of elementary 
indicators. These indicators, when converted into digital cells, are able 
to produce any sort of combinations for functional or activity maps at low 
cost; for example, "Forest Grazing Area," "Open Pit and Strip Mining," 
"Parks and Recreation Areas," or any index maps of specific environmental 
impact. This project attempts both the efficient reading of remote-sensed 
data, and the efficient manipulation of the other conventional elementary 
data. 

To bring the states along rapidly within the project calendar, the most 
gdneral compositing relationship was used — namely, arithmetic compositinj 
in which each topic map was factored by some relative weight and all were 
summed into the composite map. In some states, these composites were 
further recomposited x^ith relative weights (see Nexv Mexico and Colorado 
reports). This approach implied a knoivn relationship of the factors. The 
variety of compositing programs (GMAPS, RAP and ERGIS) was salutory and 
offered the basis for efficiency comparisons, some of xtfhich are noted in 
the state reports. A uniform data entry system was established. 

4.5 Related Activities 


A basic consideration throughout the project x^as to involve and 
communicate to as many people as possible in the region xi?ho may be potenti? 
users of the type of process demonstrated. 

As a result, a Project Briefing Paper x^as given wide distribution in 
the West, including’ copies sent to sixteen governors' science advisors. 

The project "newsletter," The Remote Sensor , was produced on eight 
occasions and distributed to the Advisory Committee and project partici- 
pants, as xvell as many other interested persons. General information x^as 
disseminated through five articles to several thousand public and private 
sector leaders in the region via the Federation Forum , a regional newslett^ 
The December 1976 Newsletter of the Rocky Mountain Science Council also 
reported on the project. 

A total of seven regional meetings were held related to project 
efforts involving some 125 persons. Another seven workshops x^ere held in 
individual states involving some 76 state and federal agency representa- 
tives. Presentations xvere given by various project participants at 
several other regional meetings and workshops and report dissemination 
and discussions xvill continue. 
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5.0 LANDSAT DATA UTILIZATION 


A principal objective of the project was to utilize digital LAND- 
SAT data for land use oriented problem-solving by combining those 
data with other multisource data via computer mapping/ compositing 
techniques. The analysis of digital data as opposed to optical products 
offers greater resolution (1.1 acres) and more efficient combination 
with other data. Colorado State University conducted all LANDSAT data 
processing, and this section summarizes their report (see Appendix A). 
Los Alamos Scientific Laboratory, with CSU, conducted the integration 
process, making the LANDSAT analysis results compatible with existing 
computer mapping techniques (see Section 6.0). 

5 . 1 LANDSAT Data 


Over 200 pages of computer listings of LANDSAT coverage over FRMS 
target areas were cataloged according to date, quality, and cloud 
cover. Calendars were then made of all possible coverages. There 
were up to eight different "image positions" which yielded data of 
interest in a given state. Each of the orbital groupings of these 
image positions had a calendar for each of the three years in question 
(1972 - 1974). On these calendars, CSU plotted all the high quality, 
less cloudy images, coding each to maintain 0, <10, <20, and <30 per- 
cent cloud distinction, as well as 0, 2, 5, or 8 quality ratings for 
each band. 

A decision was made to limit a given seasonal set of images to one 
calendar year, in order to maximize the consistency of land-use patterns 
for one season to the next. The calendar year 1974 offered the best 
choice of seasonal options in almost all cases. 1974 was also the 
closest full year choice to the timing of field selection of training 
areas in the summer and fall of 1975. 

On this basis, primary coverage selections were made for each of 
the six areas . Neither of the two common strips in Wyoming-Montana and 
New Mexi co-Colorado could be covered by a single image on any date. It 
should be noted that the coverages for Montana and Wyoming, as well as 
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Colorado and New Mexico, were interdependent; all o£ one state's test 
sites, plus one of the other's, were typically on one image, while the 
rest of the other state's sites were contained within an additional 
image. Since training data from one site were extended to use in 
classifying other sites in the same state, the dates of these two 
images had to be either the same or within 18 days^ of each other. 

Thus cloudy coverage had to be minimized over not four but eight sites 
for a given season, and cloud problems were encountered where they 
might have been avoided if the states' areas were more geographically 
distinct. 


All primary alternatives were studied, using film cassette browse 
files at CSU and USGS in Denver. IVhile the quality of the film files 
and/or viewing equipment \iras frequently so marginal as to prevent total 
insurance against unfortunate cloud locations, the search provided an 
invaluable definition of potential imagery and satisfied the objective 
that final choices were optimal. We do recommend that an upgrading of 
browse file display equipment be achieved and that a print capability 
be available at cost to the users of the regional browse files. 

For each of the following LANDSAT scenes (See Figure 5.1), computer 
compatible tapes and 1:'500,000 scale positive prints of MSS bands 5 and 
7 were obtained. 


Since the new program (see Section 5.3) for processing the multi- 
spectral tapes provided up to 18 channels of simultaneous data per 
grid cell, it was possible to use three or four calendar dates of 
imagery and, therefore, achieve a much greater discrimination of the 
land uses and cover types. Thus, seasonal LANDSAT coverage for spring, 
mid-summer, harvest season; and post-harvest season were used by four 
of the states, (Arizona, Utah, Wyoming, and Montana), and spring- svuiuner- 
harvest seasons were used by the other two (Colorado and New Mexico) . 

This posed a prodigious problem of statistical fine tuning and required 
ground truth of corresponding detail. Yet such detail was accomplished 
in only two or perhaps three of the states (Colorado, Arizona, and Utah) . 
Numerous difficulties and sources of error have been discussed by the 
CSU report and the state reports. Following are some key points: 

• First, there was needed extensive and seasonal ground truth 
for subtle differences of reflectance influenced by slope, 
aspect, moisture, vegetation overstory and understory, urban 
roofs, vegetation and mixtures in urban land uses, etc. 


^For spring coverage in Colorado and New Mexico, we were forced to 
separate dates in order to minimize clouds. Unfortunately, the two 
images did not completely overlap, and the Questa, New Mexico quad- 
rangle had its northeast corner missing. 
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OIlIGINAL PAGE 13 
PF POOR QUALmi 

FIGURE 5.1 


SUMMARY OF LANDSAT IMAGERY OBTAINED 


State 

Date 

Location* 

Quality 

Cloud % 

Test Sites Affected 




(4- Band) 


M * 

ARIZONA 

2/17/74 

39,37 

8888 

0 


.ARIZONA 

5/ IS/ 74 

39,37 

8888 

0 


ARIZONA 

8/31/74 

39,37 

8888 

10 

1 Paradise Valley, Tempe 
< (Fowler) Granite Reef 
1 Dam 

ARIZONA 

11/21/74 

39,37 

5888 

0 


COLO/NM 

S/30/74 

56,34 

8888 

0 


COLORADO 

7/05/74 

36,34 

8888 

10 

Fox Creek, Zapata Ranch 

COLO/XM 

11/26/74 

36,34 

5885 

0 


NEW MEXICO 

5/12/74 

36,35 

8888 

10 


NEK iNEXICO 

7/05/74 

56,55 

8888 

10 

Santa Fe, Quest a 

N*EK MEXICO 

11/26/74 

36,55 

8888 

0 


UTAH 

6/22/74 

41,51 

8888 

0 


UTAH 

8/15/74 

41,31 

8888 

0 


UTAH 

10/08/74 

41,31 

8888 

10 


UTAH 

6/22/74 

41,32 

8888 

0 


UTAH 

8/ IS/ 74 

41,32 

8888 

0 


UTAH 

10/08/74 

41,52 

8888 

30 


hlO/MONT 

6/01/74 

38,29 

8888 

10 


l\YO/MONT 

7/25/74 

38,29 

8888 

10 

^Decker, Hunter Mesa, 
^North-Ridge 

hTO/MONT 

9/17/74 

38,29 

8888 

0 


ItlO/MONT 

10/10/74 

38,29 

5858 

10 

Decker 

>J0\TANA 

6/01/74 

38,28 

8888 

20 


NJONTANA 

7/25/74 

38,28 

8888 

0 


^iO^TANA 

9/17/74 

38,28 

8888 

0 


MONTANA 

10/10/74 

38,58 

8888 

0 



*Path, row in BEST COVERAGE reference grid of U S. 
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• The program itself was new and complex, having been developed 
during the period when six different sets of ground train- 
ing site data were coming in replete with gaps.. 

• There arose a problem of ascertaining "retroactive" ground 
truth in the training sites, because the best LANDSAT season- 
al coverage was available for 1974, but field work had to 
begin in 1975 and in some cases go on into 1976. Only Arizona 
had near-adequate field detail for 1974 because they already 
had orthophotos for that year. 


5.2 Training Site Data 


The selection of training data from which signatures for each 
of the land-use categories are obtained is the most important part of 
any computer analysis of remote sensing data. The final classifica- 
tion maps will, of course, reflect the characteristics of the train- 
ing data used in signature development. If the training data used 
for establishing a signature for a ponderosa pine ecosystem are repre- 
sentative of a particular stand of ponderosa pine, with a specific 
canopy cover density and a specific understory characteristic, then 
the final classification maps will show very good results for this 
specific ponderosa pine type, but may not produce good results for 
other variations of ponderosa pine. On the other hand, if the train- 
ing sites for ponderosa pine were to be carelessly selected, in that 
they actually included other vegetation types, then a considerable 
confusion between these cover types could be expected. 

These problems of training site selection were magnified beyond 
their normal importance for this project, since multi date, as well 
as multispectral imagery, were used for classification. In other 
words, when imagery from several dates are combined to form signatures 
for a class, there is a tendency for the training data to represent 
only those specific fields or locations and not be representative of 
a general cover type. 

In retrospect, it is apparent that the selection of training data 
was the weakest point of this project because of the critical re- 
quirements of multidate data, and because the training site selection 
task was neither standardized nor carefully coordinated with CSU analysis 
efforts. The personnel in each of the six states had very different 
backgrounds and interests, ranging from people with very little ex- 
perience in remote sensing to individuals with a great deal of under- 
standing and competence. Under these conditions, the instructions given 
to the states and communications with the states were really not adequate 
to ensure the best results. ' 
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' ..apparent that the instructions provided the states was an 
important factor in obtaining usable training data. It is also apparent 
that these instructions, frequent and extensive as they were in most 
cases, were less than adequate and could not take the place of first- 
hand experience with this kind of problem. More thorough instructions 
are needed, perhaps- in the form of a training site selection manual, 
for future projects. Even with such instructions, the use of experienced 
personnel v-?ill always be needed to ensure the best results - '' 

T 

1 

Computerized analysis of multispectral scanner data is "dependent upon 
significant differences in the spectral reflectance characteristics of 
the categories to be classified. Basically, this means that computer 
analysis of LANDSAT MSS data is to all intents and purposes limited to 
the detection of land-cover differences. A cemetery, a city park, an 
irrigated pasture, and a golf course may all be essentially identical, 
relative to the existence of land covered by lush grassland and trees. 

The activities associated with these land uses are significantly differ- 
ent and are important to the land-use manager. A photointerpreter will 
probably have little difficulty, with large scale imagery, in telling 
the difference between the park, the cemetery, the golf course and the 
irrigated pasture. Computer analysis of MSS data, on the other hand, 
will very likely confuse these areas because of the similarity of cover 
types. 

At the initial meeting for this project (held April 8 and 9, 19-75, 
in Denver]) , the problems ivith the USGS classification system^ were dis- 
cussed, and a new selection of Level II land-use categories was prepared, 
which was subsequently refined to that shown in Figure 4.3. 

The objective was to test the absolute limits of the LANDSAT Mapping 
System. The representatives of each of the states were informed that 
crop lands, for instance, could be subdivided into specific crop types 
for consideration on this project. Therefore, each of the states pro- 
vided training data on as many subcategories as they wished. 

Suffice it to say, at this point some of the states provided many 
subdivisions (81 in Arizona) of the original 19 categories, which ul- 
timately pushed the capability of the system beyond its current limits. 
Other states kept primarily to the list developed on April 9 and did 
not make any attempt to establish subcategories. Naturally, this re- 
sulted in considerable difference in the ultimate results obtained for 
each of the states. 

The following Figure 5.2 summarizes the states' training site exer- 
cises. Every state had specific problems, some more serious than others. 
The need for standard procedures and clear communications was very appar- 
ent. There was also a need for a clear hierarchical system of classifica- 
tion for each states' land-use categories to allow for consistent class 
division or combination. 


^Anderson 1972, op. cit. 
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FIGURE 5.2 


STATE TRAINING SITE DATA SUMMARY* 



state 

Number of Land 
Use Classes 

Number of 
Training Sites 

Training Fields 
in Test Sites 

Ancillary Data 
Descriptions 

Experienced 

Personnel 

Result 

ARIZONA 

81 

245 

More 

Field Trips 
Orthophotos 

Yes 

Good 

COLORADO 

45 

161 

Yes 

Field Trips 
Existing Data 

Yes 

Good 

NEW MEXICO 

35 

48 

Yes 

Field Trips 
Existing Data 

Yes 

Fair 

MONTANA 




Field Trips 
Existing Data 

No ^ 

Fair 

UTAH 

50 

175 

i 

Yes 

Air Photos 
Field Trips 
Existing Data 

Yes 

Good 

™ming 

CCii 

19 

76 

1 

Yes 

Field Trips 
Existing Data 

No 

Fair 


*See the CSU Report, Appendix A for a detailed evaluation and inai vidua 1 
State Reports* Appendices C-H* for procedures of each state. 
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5 . 3 LANDSAT Mapping System 


Initially, this project called for the use of CSU’s pattern 
recognition system, RECOG. RECOG, however, had been developed 
primarily as an educational tool for use in the undergraduate and 
graduate remote sensing courses which are taught at CSU. This sys- 
tem did not provide for the coordinate transformations required to 
produce rectangular and scaled maps, as was desired for this project. 
Nor did this system provide for any preprocessing of the data or for 
the specialized selection and refinement of data which was believed 
necessary. As a training tool, it had been designed to provide for 
many opportunities of human interactions with the system. In other 
words, it was necessary for the operator of the program to interact 
with the data flow and make manual decisions and selections at many 
points. Although ideal for instructional purposes and desirable for 
certain kinds of basic research, it was seen primarily as an oppor- 
tunity for many human errors to occur during the processing of the 
large volume of data required on this project. 

Other capabilities not provided by RECOG, which were needed for 
this project, included means to register files of data so that multi- 
date files could be prepared for the desired multidate processing. 
Furthermore, means of displaying training data and preclassification 
analysis results were needed to permit the efficient analysis of the 
large quantities of training data received from the states. 

In summary, capabilities to be provided by the LANDSAT Mapping 
System (LMS) are listed below: 

1 

1. Fewer points of human interaction in the data processing 
flow. 

2. Improvements in the efficiency of all of the computer 
programs to be used. 

3. A coordinate transformation and scaling capability. 

4. A data preprocessing capability to include generation of 
new variables and spatial filtering. 

5. A capability to form a single registered multidate file 
for combined multivariate analyses. 

6. Improved training data analysis to include the ability 
to extract irregular fields; maintenance of the training 
data on tapes; provision for semiautomatic cleaning or 
refinement of the data; calculation of numerous train- 
data statistics; and the capability to display the train- 
ing data in map form to verify location and spatial dis- 
tribution of data performance. 

7. Provision for more versatility in classification methods 
and display of classification results. 

The following descriptions will show that most of these desired 
capabilities were met under this project. 
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The LANDSAT Mapping System is composed of four steps . Each 
performs a specific task in the production of the final results. 

The tasks involve (1) preprocessing of LANDSAT data, (2) analysis 
of class acceptability, [3) development of class signatures, and 
(4]) final classification and display. 

The descriptions in this section are based on the original de- 
sign of LMS. Some of this capability has not been achieved as yet, 
as will be discussed. A more detailed description of LMS may be 
found in Appendix A, the CSU Report. 

5.3.1 Step I 

The first module of the LANDSAT Mapping System is called 
Step I. It prepares LANDSAT data received from the EROS Data 
Center for easier and more efficient processing and modifies 
the data format to be compatible with CSU's programs. This step 
performs several data preprocessing functions, including rotation 
and rectification of the spatial coordinate, correction for geo- 
metric errors and spatial filtering. The step also provides for 
data displays in the form of conputer generated graymaps. It 
consists of four separate programs linked together by a system 
of JCL (Job Control Language) cards. The operational sequence 
is given in Figure 5.3.1. 

5.3.2 Step II 

The module of LMS called Step II was designed to make it 
. possible for 'the user to intenveave images from various dates 
and/or ancillary data to form a multidate-multivariate file of 
a particular map area. This step is composed of two programs -- 
one which combines the channels and one which trims the data 
file to a user defined map region. (See Figure 5.3.2.) 

5.3.3 Step III 

The third step in LMS is designed to analyze training data 
and compute statistical "signatures" of the classes to be mapped. 
The signatures are composed of a covariance matrix and a mean 
vector which statistically describe the spectral characteristics 
of the training fields for each class. This step is the most im- 
portant of the four that comprise LMS. The analysis strength and 
flexibility of this step gives LMS a distinct advantage over 
other scene analysis systems. It is well known that computer 
classification of multispectral data is only as good as the signa- 
tures used for the classes. The programs of Step III give the 
analyst the opportunity to optimize signature definition and 
thereby maximize the accuracy and significance of classifica- 
tion results. 

The module is composed of nine programs, each performing 
a specific task directed toward the generation of the above signa- 
tures (see Figure 5.3.3.) EXTRACT and TRANSF2 are preliminary 


- 56 - 



FIGURE 5.3.1 



IMS Step I 

r\ 

) *“4 ' 


/ 


UP TO 4 TAPES^ REPRESENTING 25 MILE 
E-W SEGMENTS OF A GIVEN LANDSAT IMAGE., 
MAY BE INPUT SIMULTANEOUSLY. 


COfjVERTS THE UNDSAT- FORMAT TAPE(s) INTO THE 
INTERNAL^ SINGLE RECOG TAPEr ONLY THE PORTION 
OF THE IMAGE NEEDED TO OVERLAY THE SELECTED 
MAP IS CONVERTED AND POOLED TOGETHER. 


ROTATE" RESAMPLES THE ORIGINAL IMAGE CELLS TO 
REPRESENT ANY SIZE RECTANGUUR OR SQUARE CELI- 
AS SELECTED BY THE USER. ADJUSTS FOR ORIGINAL 
IMAGE DISTORTIONS. SCALES IMAGE TO MAP SCALE 
(E.G.. 1:2A.OCO). 


"FILTER" THE IMAGE, 


6RAYff"L 2, 3 ... OR ALL OF THE I DIVIDUAL 
SPECTRAL BANDS IN THE ORIGINAL OR MAP OVERUY 
FORMAT. DISPUY OPTIONS INCLUDE LINEPRINTER 
AND MICROFILM GRAYMAPS. 


"LANDSAT" COMPUTER compatible tape (cCT) as SUPPLIED BY EROS 
DATA CENTER. 


"RECC6" FORMATTED TAPE (OR DISK) FILE - AS STANIDARD TAPE FORMAT 
USED THROUGHOUT THE IMAGE PROCESSING ACTIVITY. (N = 1 TO A) 
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FIGURE 5.3.2 


LMS Step II 



\ UP TO 10 RECOG FORMATTED TAPES 
/ OF A VARYING NU^eER OF SPECTRAL 
/ B,ANDS ARE INPUT. 



EACH RECOG FORMATTED TAF 
(or file) to a SELECTED NUMBER 
OF LINES AND COLUMNS DESIGNATEI 
BY THE USER> USUALLY THOSE 
NEEDED TO COVER MAP SELECTED. 
LINES AND COLUWS ARE PENUM3ERE 
BEGINNING AT LI. 


. (Df'BlNES" RECOG FORMATTED DATA FROM 
1 TO 10 SEPARATE INPID" TAPES (FILES) 
INTO 1 COWOSITE RECOG TAPE (fILE) 
REPRESENTING A MULTIDATEy MULTISPECTF 
IMAGE. 


GRAYff " I 7 2y 3. , . OR ALL OF THE 
INDIVIDUAL SPECTRAL BANDS IN CCMBINEI 
IMAGE. DISPLAY OPTIONS INCLUDE 
LINEPRINTER AND MICROFILM GRAYMAPS. 



(i + j are any integers) 
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FIGURE 5.3.3 
LMS STEP III 







. — 

mm 






GROP 



STATGE?' 


1 

cm[ - 



SIGNATURES 

/ 

/ 

Ti 

/ 



EXTRACTS THE TRAINING FIELD DATA IDENTI- 
FIED BY THE USER (RECTANGLES., IRREGULAPv 
AREAS., AND POINTS) FROM THE RECOG imCE 
FORMAT. 


HWISFORflS THE TRAINING FIELD DATA. 

FORMS RATIOS OF SPECIFIED SPECTRAL BANDS 
USES ELEVATION OVERLAYS TO ADJUST SPEC- 
TRAL BANDS FOR TERRAIN SHADO'IING^ ETC. 

"CLEANS" OUT TRAINING FIELD DATA POINTS 
V/ITH LW PROBABILITY OF BEING THE SELECTED 
MATERIAL OR HIGH PROBABILITY OF BEING 
SOTE OTHER MATERIAL, ETC. 

"GROUPS" TRAINING SETS TOGETHER HHICH WERE 
ORIGINALLY SELECTED IN EXTRACT TO RE- 
PRESENT SEPARATE MATERIALS BUT ARE NC*‘f 
DETERMINED TO BE STATISTICALLY SIMILAR. 

"STATGEN" GENERATES STATISTICS FIELDS US INC 
MAXIMUM-LIKELIHOOD APPROACH (STEPWISE 
DISCRIMINANT ANALYSIS). OTHER DECISION 
RULES CAN BE SUBSTITUTED HERE. 


"O/ERLAYS" ANY VARIABLE 
OR RESULT IN POINT FILE 
INTO A RECOG FORMAT FOR 
DISPLAY AND MAP OVERLAY. 

"SIGNATLIRES" computer 

STATISTICAL REPRESENTATK 
OF EACH MATERIAL SPECIFIF 
BY THE USER FOR USE IN 
MAPPING THE f-lATERIALS ON 
ANY DATA TAKEN FROM THE 
SAME ORIGINAL IMAGE. 

"PRIfFTS" OR "PUPES" ou 

ANY VARIABLE (s) IN TH> 
POINT FILE FOR FURTHE 
ANALYSIS IN ADDITIONAl 
PROGRAMS WRITTEN BY 
THE USER. 



POIfPT BY POINT TAPE (OR DISK) FILE. - AN INTERNAL TAPE, DISK, Aj'©/0R 
CARD FILE FORMAT miQH CONTAINS ONLY THE EXTRACTED TRAINING FIELD DATA 
AND MAINTAINS ITS CORRECT MAP OVERLAY POSITION, 












programs ii^hich respectively isolate the training data from 
their source file and expand the data to include any extra 
transformations of original variables. The resulting output is 
then ready to b_e analyzed ivith the other programs. STATGEN 
calculates various statistics on the data and is the central 
assessment tool. CLEAN and GROUP are manipulative programs 
which delete or recombine training data at the operator's com- 
mand. OVERLAY and PRIPUN are pixel-by-pixel output devices, 
the first in map form and the second either print or punched 
cards. SIGNIT generates the class signatures when the training 
data has reached a satisfactory level of performance. 

5.3.4 Step IV 


The final step in the LMS system performs the classifica- 
tion of each pixel of an image area and then displays that area 
via some visual medium. At this time the only visual media that 
have been developed for this step are microfilm or a printed sym- 
bol map. This step uses the signatures produced by SIGNIT in 
Step III, and the target data file produced by Step II to generate 
the desired output classification map. 

5.3.5 Present Status of LMS Development 


The LANDSAT Mapping System is a valuable system because 
of Its many capabilities and options; but, because of the problem 
discussed below, the system is in many ways difficult to use and 
in some cases more expensive than necessary. A complete docu- 
mentation of the system is not available due to time and money 
constraints and is very badly needed. 

During the development of LMS, the programs were written 
by a small staff under very tight time constraints. Therefore, 
only minimal documentation was written. Also, because of the 
lack of time, parts of the system are not efficient. A particu- 
lar case in point is program OVERLAY in Step HI, which in most 
cases is too expensive to run, but remains a valuable tool in 
determining the character of individual results during Step III 
processing. The package has not been entirely completed (90 
percent complete) and contains some problems which currently 
impede its use. 

5.3.6 Future L^fS Development 

The LMS package, because of its many options and capabilities 
in perfonning classification of various areas, is a valuable tool 
in the field of remote sensing. But as it currently exists, it 
would be difficult for anyone to use the software because of the 
lack of optimization and documentation. 
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The decision points encountered in the operation of LMS 
need clarification in terms of the implications of alternatives 
presented. In order to export LMS, subsequent to the above im- 
provements, the preparation of manuals for training data collection, 
software operation, installation, and program maintenance must be 
achieved. These have been outlined in Appendix A, the CSU Report. 


5.4 Processing the Data 


5.4.1 Preprocessing LANDSAT Data 

In this section ive will describe the problems encountered 
and the results obtained while using Steps I and II of LMS for the 
preprocessing of LANDSAT data. Unfortunately, Steps I and II of 
LMS were not fully developed when we initiated data preprocessing 
for several of the states. Therefore, we had a continually changing 
system ^^rith little or no documentation during this preprocessing stage. 
This compounded and greatly magnified the normal problems always en- 
countered when using digital data. First, we will discuss the kinds 
of problems that any user of LMS will likely encounter, and then we 
will briefly describe the problems and results associated with the pre- 
processing of data for each state. 

There were many problems which resulted from use of software 
which was under development and undocumented. They were very frus- 
trating and caused many delays and some cost inefficiencies. These 
problems should not be encountered by future users of LMS and so 
will not be discussed here to any extent. The greatest importance of 
these problems lies in the impact they had on the design of “Steps I 
and II of LMS. A very regimented step-by-step procedure was required 
if this system were to be used by relatively inexperienced personnel. 
Initially, a series of equations were written for determining the 
location of data for target areas on the digital tapes. The use of 
these equations resulted in many errors and many repeat operations of 
Step I. Step I was modified, therefore, so that the user need only 
identify one corner of the area to be extracted. The software was 
then designed to compute the other parameters required to extract the 
correct data for the target area. This greatly reduced extraction 
errors but did not completely eliminate them. 

Mistakes in extracting data from the digital tapes continues to 
occur, primarily because of a lack of exact correspondence between the 
region shown on the 1:500,000 scale printed images and the region en- 
compassed by the data on the digital tapes. In other words, there were 
instances where the image prints did not correspond exactly with the 
data on the computer compatible tapes. These vertical location errors 
were at times as great as 50 lines. Horizontal location errors were 
generally of small magnitude and did not present a serious problem. 

These difficulties could be circumvented, of course, by extracting 
data for regions which extend significantly beyond the target area. 

This can be costly since more data is processed than necessary; but 
in the final analysis, it will probably save money since it will elimin- 
ate the need for repeat computer runs when target areas are missed. 
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Radiometric errors, caused by differences in the sensitivity 
and calibration of detectors on LANDSAT, produce marked, horizontal 
striations on the images. This problem is particularly serious for 
MSS band 7. The rotation and filtering of the data reduce the serious- 
ness of these radiometric errors, but for some of the ^^^orst images, the 
problem was still quite significant and can be seen on the final classi- 
fication results. The LANDSAT tapes do provide some calibration informa 
tion which could be used to further reduce radiometric problems. At 
present, IMS does not provide an algorithm for making use of this 
calibration data, so no additional correction i?as possible under this 
project. Provision for such correction will be a high priority item 
for continued development of LMS. 

Parity errors and the associated missing data affected the images 
for several states. It was not a serious problem except for one Wyoming 
quadrangle which will be discussed later. When only one line of data 
is missing, the two-dimensional filtering of the data provides some fill 
in capability. This is far from adequate for the filtering intensities 
used on this project, however, and a surface fitting algorithm should 
be developed to replace missing data. 

Unfortunately, it is very difficult — if not impossible — to 
identify radiometric errors and parity errors on the image prints used 
to select image dates. For this reason, the quality of the data is 
unknown until preprocessing begins. 

Another data problem v^hich was not detectable on the image prints 
was the presence of thin cirrus clouds or small cumulus clouds in some 
of the images. Some very thin clouds in Wyoming were not detected until 
their effect was noted during the analysis of training field data. It 
would be desirable in the future to develop software to adjust radiance 
values under thin cirrus clouds, based on the radiance of near-by areas 
and to fill in data for locations underneath very small cumulus clouds. 

The same filter weights (filtering intensity) were used for all 
image dates for all states. This was undoubtedly not optimum filtering. 
In those areas, such as urban scenes, which contained many small but 
important features having high spectral contrasts, it would have been 
preferable to filter with less intensity. For certain wildland scenes 
and some large water body areas, however, it would probably have been 
desirable to filter with greater intensity. The intensity used was 
probably adequate for most situations with the possible exception of 
the urban, suburban areas. 

No problems were encountered in using the filter algorithm, but 
it is likely that improved results could have been obtained in some 
quadrangles if the algorithm had been more flexible (provided different 
weights for different scene conditions) . 
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The between-date registration of scenes and the overall 
registration of the scenes to the 7.5-minute quadrangle maps presented 
very few problems in agricultural and urban, suburban areas., but 
proved to be quite difficult in most wildland regions. The serious 
problems encountered generally resulted from either changing scene 
characteristics from date to date, or the lack of image contrast for 
certain quadrangles. 

For example, the Dromedary Peak quadrangle- was by far the most 
difficult Utah quadrangle to register, but it was not because of lack 
of scene contrast. The problem lay in the changing characteristics 
of the scene in this very nigged terrain. Specifically, the apparent 
location of ridge lines changed from date to date because of the shift- 
ing of shadow lines as the sun zenith angle changed. 

On the other hand, in Wyoming and Montana some of the quadrangles 
contained no major topographic or geographic feature for use in estab- 
lishing registration. These relatively flat areas contained no ridge 
lines, no lakes, nor any other feature to ensure accurate registration. 
Under these conditions, there was probably a certain amount of residual 
errors in the final registration results. In most instances, however, 
a careful examination of the imagery located several vegetation or 
topographic features which were consistent over the dates in use. It 
is likely, therefore, that the between-dates registration was, in most 
instances, within one pixel. Final registration to the topographic 
map, however, could have contained errors of two or three pixels for 
several quadrangles. 

5.4.2 Analysis of Training Data 


Step III is that portion of LMS which provides analytical tools 
for the evaluation and optimization of training data. This could be 
called the signature development step of LMS. Without question, this 
is the most important and critical step of any image analysis system, 
since the final accuracy and quality of the image analysis are dependent 
upon the accuracy and validity of the signatures. The computer analysis 
of the image data can do nothing more than compare the spectral charac- 
teristics of each pixel with the spectral signatures which have been 
established in Step III. If the spectral signatures are wrong, the 
classification results will also be \vrong. 

Since the training data generated by each date differed appre- 
ciably in character, quality, and documentation, treatment of those 
data varied accordingly. 

Summary of Step III procedures: 

1. Extraction of Training Data . The first procedure in Step III 
calls for the extraction of training data from the master 
files created in Step II. This is accomplished with program 
EXTRACT whose output is in the form of a POINT file where each 
individual pixel within the training fields can be identified 
and analyzed. 
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2. Generation of Nei"? Variables . After extracting the 
original data from the master files, there is the 
option of generating new variables through transforma- 
tions performed on selected bands from the original 
data. This is accomplished with program TRANSFORM. 

For this project, new ratio variables were generated 
by ratioing MSS bands 7 and 5 and bands 5 and 4. - 

3. Initial Class Grouping . The user now has the option 
of reordering the training data pertaining to speci- 
fic classes to facilitate the analysis of a subset 
of these classes. This is accomplished with program 
GROUP. This grouping of classes is needed because 
program STATGEN analyzes a user specified quantity of 
data, always beginning at the start of a point file. 

If only part of the classes is to be analyzed at any 
given time, therefore, it is necessary that these classes 
be at the front of the POINT file. The reasons for 
analyzing a subset of classes are discussed later in 
this section. 

4. Analysis of Training Data . Program STATGEN generates all 
of the statistics used in analyzing the training data. 

These statistics include the mean vectors and coverian'ce 
matrices of classes, a measure of the value of each of 

the variables in terms of F ratios, an F matrix for measur- 
ing the separability of each class from every other class, 
a within-groups correlation matrix to identify the correla- 
tion between variables, and a list of a posteriori proba- 
bilities used in the nest procedure to eliminate 
individual pixels which are degrading the reflectance 
characteristics of specific classes. 

(At this point, many options are available to the user. If all the 
statistical measures of class separability are satisfied, he can return to 
Procedure 3 above and access and analyze the next subset of classes. If 
the class performances are substandard, he may choose one or both of the 
next two following procedures:) 

_5. Cleaning Training Data . Program CLEAN is used as noted 
above to delete anomalous pixels in selected classes . 

In many instances, it may be .desirable to operate pro- 
cedures 4 and 5 in an iterative fashion, until the per- 
formance of each of the classes is acceptable. In other 
words, as pixels are deleted, the mean vectors and covariance 
matrices (this is the signature for each class) are modi- 
fied and the classification results are usually improved. 

6. Combining or Deletion of Classes . At this point, program 
GROUP can be used to combine classes which are not separable 
to establish a lower-level class or at least a class that 
is a combination of previous individual classes. This pro- 
gram can also be used to delete any class from subsequent 
analysis. 
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7. Final STATGBK Run , The final step in the detailed analysis 
of the training data normally will call for a combination 
of all classes used in each of the subsets of classes for 

a combined computation of classification statistics. This 
is to ensure a minimum of confusion between class subsets. 

8. Calculation of Signatures . After all of the analysis and 
manipulations performed have been concluded to the satisfac- 
tion of the user, program SIGNIT is used to create final 
signatures for each class for use in Step IV. 


5.4.3 Processing of LANDSAT Data 

Upon the completion of the development of signatures for each land 
use category, processing of the LANDSAT data for each quadrangle to be 
classified began. This was accomplished in Step IV of LMS, using programs 
MAP and DISPLAY. The inputs to program MAP are (1) the mean vectors and 
covariance matrices which define the signatures for each of the classes^ 
(2) the matrix of multispectral and multitemporal data for the area to be 
classified. Since a normal distribution of spectral values is assumed, 
program MAP calculates the probability density function for each pixel 
using the equation: 




C2 H) 


where }C is the matrix of data in the spectral, 
temporal and spatial domains; 

N is the san 5 )le size;. 

li- is the mean vector for class C- ; and 

Z-j is the covariance matrix for class Cj|^. 

This equation is calculated for each class having a finite probability of 
occurrence within the quadrangle being mapped. For each pixel a decision 
for one class versus another is based on the simple comparison of density 
values given by 

^ I ^i ^ s 1 (decide C^) 
f (X^ I C 2 ) 

There is always the probability, of course, that a given pixel may 
not belong to any of the categories or classes being considered. A 
probability density threshold T is used, therefore, to establish a lower 
threshold, below which a decision will be made to leave the pixel 
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unclassified. This decision algorithm may be expressed by: 

If f (Xj I C^) < T (decide none) (5.3) 

for all i 

The threshold values T are selected on the basis that the probability 
density values for normally distributed populations of variables will 
have a Chi Square distribution. Therefore, the values for T are selected 
from a Chi Square distribution table. (See Figures 5.4.1 and 5.4.2 for 
examples of the 0% and 10% threshold land use classification maps.) 

Under normal conditions, the processing of LANDSAT data (after a 
proper selection of classes and an adequate development of signatures) 
becomes routine. Such was not the case for this project. As already 
noted in previous sections, the selection of classes and the selection 
of data for the development of signatures left much to be desired in 
all of the states. In addition to the problems generated by these short- 
comings, two new difficulties arose to compound the problems. 

First, the covariance matrices would not invert, as is required 
by equation (5.1). Second, there were money problems since the large 
number of classes and variables raised the cost of classification beyond 
budget limitations. These problems are briefly discussed below, along 
with the procedures used in processing the data. 

Frequently large matrices will not invert because of one or the 
other of two common problems. First, if a row or column of a matrix is 
repeated, the matrix will not invert because of the complete correla- 
tion between the identical columns or rows. Secondly, a large matrix 
will often be noninvertable when there is a large difference in magni- 
tude of data values within the matrix. Essentially, this second problem 
arises from a lack of sufficient accuracy in computed values, resulting 
from limitations in computer word size. 

Both the problems noted above appeared to exist for the signature 
covariance matrices developed under this project. It was discovered, 
for instance, that none of the matrices would invert when ratio variables, 
generated by taking the ratio of original LANDSAT variables, were employed. 
This appeared to be due to the correlation between the generated new 
variables and the original variables from which they were formed. On 
other projects we had successfully used ratio variables, but these projects 
used imagery from only one date, thereby reducing the total number of 
variables and the size of the matrix. Furthermore, the previous projects 
only involved classes of natural vegetation types. A comparison of the 
correlation between ratio and original variables for natural vegetation 
classes and such classes as residential and agricultural crops has shown 
a higher correlation for the latter classes. 
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FIGURE 5.4.1 
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FIGURE 5.4.2 

Partial Map: LANDSAT LAND USE CLASSIFICATION -- 10% Threshold 
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SCALE 1:24,000 


LANOSAT IMAGE OATES: Miy 15, 1974. RowemUcr 20. 1974 
duty 7. 1974 


CLASSIfICATION SUH 


(Cla»>tfiCAtton Thf«tHtit4 10\) 


CLASS DESCRIPTION 

PIMVON'JUNIPLK • open 
9Sf/SJ tuth unJemvAth 
ftNYON -JUNIPER • open 
90P/10J toll iMMlAmrAth 
flNYaN-JUNiPEJI • open 
70P/30J soil ^ndtrnoAth 
riNYON-JUNIPER • clostd 
covoA'very Uttto toil 
PINYON-SACE • SO AA|o/ 
to plnyon 

HUCO CONIFER • »p«n 
lOP/50J/20PP/700i|F 
MIXED CONIFEP - closoJ 
2SJ/SOPP/2tOiF 

corroNMooo 

MtSTtRN MlEATCRASt 
H1.A00W 

GOOD CONOITIO*! SAGE 
PAIR CONDITION SAGE 
POOR CONDITION SAGE 
NINTCRFAT 
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TAICING PONO 
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un 
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240 

704.1 

SOI 
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I9S 

270 


CUSOIFICD ZiZlfi 

UNCUSSIMLD ISIOO 
total ARU 3IJ7S 


I. THESE LAND USE CUSSIFICATIOnS NTRl HAUC IV 
CUMPUrCR PkOCtSSUlC OF lAMiSAT OCW-IMY . 
MUCTlUATi: FlCiS MRl PkOCLSsCO AS A si>a:il 
MJCTIVARUT li tlAJE. 

2. THE DATA MINE FILTl. J USING A TH) OINLNSKKAL 

NEigrrED averaging I.CIMIGUE this hit fkOUUCt 
tKKOkS OK UNCUSSirilU BUKDEHS BIT NUN CLASSU 
HAVING URGE RAIUOU IKIC OlFF'LHrNClS. 

1. CAOI SmiOL kEPRlsmtS AH ARIA OF 1.1573 ACKI.S. 
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DISIMlIiUnuN. TIU TUlintlOlO PiRCIS'lAi.C MHIID 
It A TRUt PKObAlUirY THRIUMD. 

Tilts MAP PRirAREO $Y: lOLORAIO SIATE USIVIR'IU. 

FORT COII.IRS, CUU7RA00 AOS^S. MM>IKOOMUAi | 

TO IIHIKAIUIN M Atl i) Ik.tl'.IAIK SIAM I'|..),H 
NiiMitrr NAb s*7.'l)l. Fui AJJiIiuujI llll..nkAlltMi 
Coiilact: Lufriit L. M4*ii»||, Dlrai:i Jr*RcM>*la 

ScnAlny liogrdmi. (303) 491 t|47 or 4V|.Si.bl 


Whatever the ultimate cause, it was determined that the matrices 
would not invert for some classes when ratio variables were included. 

All such variables were, therefore, eliminated for the processing of 
the LANDSAT data. This undoubtedly reduced the separability of certain 
natural vegetation classes, since the separation of almost all rangeland 
and certain forest type classes are improved when ratio variables are 
employed. ^ 

For certain classes, we also found that covariance matrix values 
covered a range from less than 10 to several hundred. This occurred 
because during some seasons there would be very little variability in 
certain spectral bands, whereas for another time of the year^ the radi- 
ance values might become highly variable. 

These problems were of a very serious nature because whenever 
the covariance matrix failed to invert properly, the classification 
results were completely meaningless. This being the first project for 
which we had undertaken multitemporal analysis of data, we had not pre- 
viously encountered such difficulties. Not anticipating these problems, 
we had not allowed time for their solution; furthermore, since other 
difficulties had already caused project delays, we found ourselves with 
very little time to seek an optimum solution. Within the time available, 
the only solution was to undertake a modification of the covari- 
ance matrices. The first modification has already been noted; this was 
the elimination of all ratio variables. The second modification essen- 
tially involved the modification of signatures. 

From the analysis of the training data, it was determined that 
the sample size and the lack of spatial distribution of our samples 
were probably, at least in part, the cause of both extremely small 
covariance and extremely large ones. This is merely a matter of having 
a sample size too small to adequately estimate covariance values. This 
statement is supported by the fact that these problems did not occur in 
those states and for those classes where a very large, widely distri- 
buted set of training data was available. 

When the problem did occur, it was solved by computing a common 
covariance matrix for a group of classes having a similar reflectance 
characteristic. For instance, all of the forest type classes were 
grouped together (when necessary) to form a common covariance matrix. 

In this way, the sample sizes were increased and an invertible covariance 
matrix was obtained. 

Better results would undoubtedly have been obtained if there 
had been time to establish larger sample sizes and a proper selection of 
classes themselves, in seeking a solution to these problems. 

The cost problem associated with too many classes, and too many 
variables required some extraordinary modifications of our classifica- 
tion algorithm. It has been noted that equation (5.3) can be used to 
set a threshold for discriminant values, below which no class will be 


^laxwell, E.L. Multivariate Systems Analysis of Multispectral Imagery** 
Photo. Engr. and Rem. Sens .. Vol. XLII, No. 9, pp 1173-86, 1976. 
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selected. Using the same assumption, that the discriminant values 
should have a Chi Square distribution, the decision function given by 
equation (5.2) was preceded by a decision function 


£ (X^ j C^) > T* (decide i) (5.4) 

where the threshold values T' were selected for a 75 percent probablity 
that the discriminant value belonged to the given Chi Square distribu- 
tion. Using this combination of equation (5.4) and (5.2) whenever a 
class satisfied equation (5.4), the decision was made at that point and 
no additional classes were considered. If no class was sufficiently 
strong to meet the requirements of (5.4), all classes were analyzed and 
the class decision was based upon (5.2). 

This decision algorithm was effective in reducing the cost of 
classification only if the most probable classes for a given quadrangle 
were analyzed first. Therefore, small samples of data were used to 
determine the dominant classes for a quadrangle. Specifically, pixels 
from every ninth column and every ninth line were selected for the 
initial analysis to determine the class order to be used in the final 
data processing. This is called a 9x9 test. After this determination 
was made, we also performed 4x4 tests (using every fourth column and 
every fourth line) to determine which classes could be eliminated from 
the final processing for a given quadrangle. This test was based on 
the hypothesis that if a class were not selected for a six percent san^le 
of the quadrangle data (evenly distributed over the quad) , it was un- 
likely that class had a significant population within the quadrangle. 

In effect, an a priori probability for the occurrence of classes was 
being determined. Those classes whose a priori probability was zero 
(based on a 4x4 test) were eliminated from the final processing of the 
quadrangle. 

After all of the tests had been performed and the best selection 
of classes and class orders had been established, the final processing 
for each of three quadrangles within each state was undertaken. In some 
instances, the thresholding decision function, given by equation (5.4), 
was not used when it was determined that its performance was giving 
inconsistent results. 

5.5 Results of LANDSAT Mapping Verification 

A detailed verification analysis and comments about the classifi- 
cation results for each state are, of course, contained in the final 
reports prepared by each of the states (Appendices C through H) . This 
section will comment briefly on procedures, the apparent classification 
results based on training data, problems encountered and obvious 
shortcomings of the results. 

To obtain an accurate assessment of the quality and accuracy of 
any data collection or information gathering effort requires carefully 
designed and statistically valid sampling and analysis methodologies. 

Since each state undertook its own verification program, there exists 
a wide variation in the methods used and the accuracy of the data 
collection used for verification. 
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Ideally, a random selection of sites, from which verification data 
is obtained, should be used. Any other method of data selection will be 
found to produce a bias towards certain classes or certain field .con- 
ditions. Unfortunate 1-y, the time and money required for such an ambitious 
assessment was not available. CSU prepared an overlay map showing 10 
acre cells, which had been randomly selected by a computer program, for 
use by each of the states. Even under the best of conditions, it was 
not possible to obtain a sample size as large as would be desired. 

However, it was realized that the random sampling scheme would not 
be applicable to certain land cover categories such as small area urban 
or industrial land uses. For these categories the best approach was to 
identify and verify most of the classes. 

Methods used to collect the ground truth for verification varied 
widely, from vegetation type maps of the Forest Service, to photointer- 
pretation of aerial photography, to actual onsite visits. Because of 
this variety of sources and the associated problems of accurate geo- 
graphic referencing (for the 1.1 acre cells), the quality of verifica- 
tion data was variable. In some states the quality of the verification 
data was collected at the beginning of the project, when training data 
was collected, and in other states the verification data was collected 
after the computer classification mapping was finished. Neither method 
is entirely satisfactory: a simultaneous collection of verification 

and training data probably biases the verification results toward a 
favorable report; whereas collection of the data after computer classi- 
fication mapping is completed may result in a more negative result than 
is deserved, because of differences in criteria for selection of class 
fields, and the time elapsed since the imagery date. 

Accurate ground location of individual 1.1 acre LANDSAT pixels 
also presents a constant problem of verification accuracy. It would 
require a very time-consximing and expensive survey of verification 
site locations to know for sure that the mapped pixels and the veri- 
fication spots are in the same location. This is one reason 10- acre 
cells were used for verification procedure. 

In order to obtain a more complete analysis of accuracy, it would 
be desirable to assess the accuracy, both at the highest level of 
classification and at loiter levels, as classes are aggregated. 

States did not assign statistical confidence limits to their 
verification results. It can be said that, almost without exception, 
sample sizes were too small to achieve a high degree of confidence in 
verification accuracy. Basically, the difference of sampling schemes 
and analysis methods made it impossible to consolidate results from 
state to state. In reading the state results, as summarized in Figure 
5.5, careful consideration should be given to the sample size. The 
technical literature indicates that a sample size of at least 50 is 
desirable. Since, in many instances, the accuracy figures are based 
on sample sizes of less than 10, we treat the results with a great 
deal of caution. Nevertheless, some general indication of the per- 
formance of the LANDSAT analyses is provided. 
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FIGURE 5.5 


STATE VERIFICATION PROCEDURES AND RESULTS 


STATL VUUI ICATION PllOCI DUUfS AND R1,SUI,1{> 


State 

Simple 01 Complex 
C.itejjoiiCh lJi>oU 

Verification Sample Distribu- 
tion Method (Stilt I St j cal ly 
Random or Convenient Use) 

Intel pi et. It ton and 
8LalisLic.il Methods 

lx.iiiipie Rela- 
tive lljgji, 

Sloj 1 ng 
fat ego i ICS 



ARIZONA 

Veiy many (81) 
catciioiioh, yet 
some are com- 
plex. Selection 
Ubcd also season- 
al changes. 

Used arbitiaiy sites based 
on much iivailablo field 
data. 

Lxlensjve tieatmeni. 
Used diagonal luatiix, 
detaiJcd discussion. 

Li ti us 
Ol chard s 
Vi ncyaids 

(h il])CS 

Residential 
Pastille 
Mes(|in le 

See Section 3 2 ot Auz 
report, Appendix C, for 
extended statistical and 
procedural coverage. 

COLORADO 

Well stratified 
forest, brush § 
meadowland cate- 
gories. 

Began with random over- 
lay selected about half 
of sites, enlarged some 
sites 

Set up computerized 
verification procedure 
for tliree types of 
eirors Pioduced di- 

agonal matrix 

Meadow 

Coniferous 

Deciduous 

Range 

White fir 
Douglas Tir 
Certain shrubs 
§ grass 
(Irrigation) 
Ponderosa/ 
Meadow com- 
bination 

See Section 3.2 and Appen- 
dix 8 of Colorado Report, 
Appendix D for extended 
statistics and evaulation 

MONTANA 

Stratified 

categories 

Used random site over- 
lay and extensive use 
of "V'' forms, aerial 
photos 

Simple frequency tabic, 
with narrative evalua- 
tion of probable causes 
of error 

Diy hay 
Dry cereals 

Mines 

Wheat 

Pontlorosa hay 
Ii u gated 

See Section 3 2 of Mon- 
tana Repoit, Appendix C. 

NEW MEXICO 

Many categories 
used — some 
complex, some 
functional. 

Used given random sites 
overlay, "V*’ forms, 
field visits 

i 

No statistics reported 

Sagebrush 

Meadows 

Confifers 

Irrig. wheat 
Low Dens. Res. 
Commercial 

See Section 3 2 of New 
Mexico Report, Appendix V 

UFAM 

Well stratified 

Data collected early along 
with training data collec- 
tion, used three methods 
of site selection. (1) Ocu- 
lar examination (LAND SAT 
map with conventional maps 
Q photos), (2) pre-selected 
sites, (3) random sites. 

Used aerial photos ' 

Used diagonal erior 
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These verification indications cover many procedural steps, 
from initial class definitions to training site selection, recording 
of land cover variables on a retrospecjiive basis (in aH states except 
Arizona) , development of the new computer mapping program, handling 
multidate imagery, and numerous variables inherent in verification site 
selection and description. Indeed, verification procedure should be 
treated as a study in itself, in order to specify what is needed for 
an operational regional LANDSAT updated data bank. The entire subject 
of ground truth needs vigorous control. The six*-state experience with 
this problem has led to certain recommendations by the Federation, 
particularly for the establishment of a system of Standard Test Sites . 

5.6 State Classification Synopsis 

Because individual states were given the option to test a variety 
of classif ication and verification procedures, the results are not 
consistent from state to state. 

5.6.1 Classification Results From Arizona 


On the basis of Level I detail, the results for Arizona 
appeared reasonable and consistent. The aggregation and simplif lo- 
cation of desert shrub categories resulted in more informative 
and useful results than would have been the case if a finer 
category structure had been chosen at a cost of reduced accuracy. 

It was apparent that rock classes were being confused with 
such development classes and shopping center, mobile homes, and 
the higher density residential categories. This isnot surprising 
since most of these represent non vegetated land. 

Agricultural classes were an exception to the maintenance 
of areal consistency. Many fields are shown with unlikely com-- 
binations of several crops. Field boundaries are usually inaccu- 
rate for two reasons. First, the problem of unavoidable within- 
pixel mixtures resulted in misclassif ications . Secondly, spatial 
filtering tends to spread the boundaries between areas having 
different radiance characteristics. 

Among the residential classes, the single family unit (SFU) 
types showed surprisingly appropriate distribution with the 
exception of commission error into barren areas. Commercial and 
shopping center classes performed in a similar manner to SFU’s. 
Mobile homes were generally classed correctly wherever they existed, 
but the class was seriously overextended into other areas, pro- 
ducing large commission errors. Power plants and oil refineries 
were enigmas throughout the analysis, and showed expected incon- 
sistency. Airport was, for practical purposes, entirely replaced 
by mobile homes. The factory class did not perform well at all, 
appearing in anomalous localities such as highway interchanges. 
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5.6.2 Classification Results From Colorado 


Since CSU conducted the verification program for the Fox 
Creek quadrangle in Colorado, more information about the procedures 
used (their strengths and weaknesses) was obtained than for any 
other state. 

The agricultural categories were distributed erratically 
due to the same problems of multidate definition which *« 5 ere 
observed in Arizona — with three dates of measurements "^defining 
classes whose phenologies were at best distinct on only one date, 
the classifier tended to locate the categories according to com- 
pound field conditions rather than actual crop type. There were 
also uncertainties arising from questionable field descriptions, 
because of dependence on farmers^ memories for ground truth. Wet 
pasture was distributed too widely being confused with small 
towns and possihley other crop types, which were omitted from the 
training program. 

The valley floor natural cover classes seemed to perform 
quite well. Basalt butte, alkaline bare soil, greasewood, 
riparian grass, and cottonwood all were located in coherent and 
reasonabley distributed patterns. 

Development categories were successfully mapped largely 
to the extent that training areas comprised a large portion of 
their distribution. In other situations, small towns were found to 
replace gravelly areas on alluvial terraces, while the town of 
Manassa was almost entirely mis classified as pasture, alfalfa and 
plowed field. In addition, commercial classes were cropping 
up along roadways, possibley due to turn-out widenings. In any 
case, outside of the training areas, our development classes 
were very unpredictable. 

In summary, the substandard results obtained for agricultural 
types could be improved through better multiseas onal ground truth. 
Development categories were poorly devised, and reselection by 
strictly visual criteria could improve results. In the natural 
cover categories, encouraging results were obtained for a majority 
of the classes and high error rates in some of these natural classes 
were caused more by faulty category selection than any other system 
problem. 

5.6.3 Classification Results From Montana 


By using a combination of supervised and unsupervised se- 
lection of training data for Montana, reasonable classification 
results were obtained* They undoubtedly do not represent the best 
results that could be obtained from LANDSAT multidate processing, 
however, and this fact should be kept in mind when reading the 
Montana state report. 
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5.6.4. Classification Results From Mexico 


The Nex7 Mexico mapping performance ranged from satisfactory 
to poor, largely controlled by the degree to which training signa- 
tures were extended geographically. The Santa Fe quad in particular 
was substandard, due to the fact that all of the training signatures 
for natural cover came from fields approximately 50 miles to the 
north. This in turn was the source of incompatibility between 
the extended natural types and the locally trained development 
classes, resulting in persistent mislocation of development 
classes on natural areas in the Santa Fe quad. The opposite 
extension problem was apparent in the Taos quadrangle. Old and 
new medium density residential seemed to be erroneously mapped 
over gravelly alluvial surfaces, and on forest-range boundaries. 

Agriculture in New Mexico was represented by either irrigated 
meadow, alfalfa, or irrigated wheat. The latter was under- 
represented, and confusion was high betx<reen meadow and alfalfa. 

In the northern quadrangles, the predominance of natural cover 
classified by local signatures yielded results which were encouraging 
at the Level I/II detail. Pinyon vs. grass/shrub range was particu- 
larly well differentiated according to comparisons of results vs. 
forest tinting (30%) on topographic maps. Within the forest types, 
major species breakdowns seemed to distribute themselves properly 
along altitudinal and aspect variations. Herbaceous classes were 
mapped with noticeable correlation to that portion of the moisture 
gradient which their training fields exemplified. 


Guadalupe Mtn. was the highest map of the three, with the 
only serious problem arising from a lack of training on the Rio 
Grande gorge. There is some question as to exactly what information 
the sage and piney on- juniper breakdowns convey, but they nevertheless 
differentiate between real ground conditions, and to that extent are 
inherently valuable. 

In summary, the mapping performance within a class was so 
variable from quadrangle to quadrangle, that any single conclusion 
or accuracy estimate would have to be qualified by regional 
variation descriptions. Verification of the Nex^ Mexico results 
will be complicated manyfold by this need for careful stratification. 

5.6.5 Classification Results From Utah ^ 

The classification results for Utah appeared to be quite 
reasonable and satisfactory. It x^ras obvious that in the urban 
areas the residential, commercial, and industrial classes were 
being confused among one another and with certain other classes. 

One of the industrial classes, for instance, showed up in the mud 
flats of Salt Lake, where no industry existed. Apparently, the 
reflective characteristics of the mud flats was similar to con- 
crete and macadam. In the final analysis, some urban classes were 
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not used in processing the data since it was apparent that the 
accuracy of classification was going to he too low to be of any 
value * 


Classification of agricultural crops and natural vegetation 
cover types appeared to be xn good order. In the Dromedary Peak 
quadrangle, for instance, it was noted that the different vegeta- 
tion types seemed to be found primarily in fairly large homogeneous 
regions . 

The division of the classes xn the Dromedary quadrangle into 
subclasses based on slope and aspect was quite successful. By assign- 
ing the same symbol to the subclasses , the final map products looked 
quite reasonable, based on the information obtained from the train- 
ing data. 

CSU found it necessary to use modified sxgnatures (average 
covarxance matrices) to obtain satisfactory results. 

5.6.6 Classif Icatxon Results From Wyoming 

The initxal classifications for Wyoming were very poor. It 
was apparent even from a preliminary examination that the selection 
of classes and the selection of training data had not been adequate 
to achxeve satisfactory results. 

It should be noted that the Wyoming report was based on these 
initial results and does not, therefore, represent an assessment of 
the true potential of multidate LANDSAT processing. This was parti- 
cularly true since the key quadrangle for Wyoming, the Buffalo 
quadrangle, was probably the poorest of the three that were 
processed. 

Because of the obviously bad results obtained for the 
Buffalo quadrangle, a decision x^as made to obtain new traxnxng 
data for a nex7 selection of classes for that quadrangle and 
to reprocess the data. This reprocessing was completed just a 
short time before the completion of this report and shows very 
satisfactory classif icatxon results. Unfortunately, these re- 
sults were not available in time to be included in the Wyoming 
state report. 

The classes for the repeat analysxs of the Buffalo quad- 
rangle xncluded nonirrigated and irrigated rangeland and hay 
fxelds segregated according to productivxty or biomass values. 

These production classes or range condition classes performed 
quxte well and show the effects of over-grazing and availability 
of xrrigation water. 

All of the classes xv^ere xncluded in the final processing 
of the Buffalo quadrangle; no signature modifxcatxon was necessary 
and only equation (5.2) was used for a decision algorithm. 
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6.0 COMPUTER COMPOSITE MAPPING 


6.1 Approach 

A complementary relationship between LANDSAT maps and other 
.sources and forms of land use information is made possible through 
grid mapping and compositing. Such relationships for analysis and 
simulation are particularly useful if not essential, for full 
utility of these two' powerful technologies in regional planning, in 
resource studies for agriculture and mineral industries, etc. 

This project used several procedures of cellular mapping to 
siq)plement LANDSAT CCT products. It was necessary to map from other 
sources certain land use ''activities" which could not effectively be 
mapped by LANDSAT and had no distinguishing reflectance character- 
istics — essentially only man-made patterns. For a spatial data bank 
there is available a wealth of data such as population and housing 
census, building and zoning maps, agricultural crop and soil maps, 
utility system maps, etc. Most of this data is convertible into 
whatever grid cell pattern is adapted for LANDSAT CCT maps. In such 
flexible form, it is possible to reconstitute some of the complex re- 
lationships existing on the ground and in social and economic activity, 
which has been "abstracted" in various ways and through various media. 

6.1.1 Compositing 


IVhile data banks require many sources and forms of infor- 
mation, the optimum use of LANDSAT CCT maps in such a bank is to 
bring in the most current synoptic land surface state. For example, 
the several exercises in "urban construction feasibility mapping" 
carried out in Arizona, New Mexico, Utah, and Wyoming place most 
weight on maps of soils, slopes, and flood-prone areas, relying on • 
LANDSAT mainly for constraints such as prime agricultural lands, 
drainage ways, etc. to be protected from urbanization. (See Figure 
6 . 1 . 1 ,) 


In general, computer mapping and manipulation processes can 
be divided into two major categories: (1) the grid or cellular 

system, and (2) the polygon system. ' The cellular system can be 
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FIGURE 6.1.1 


USING A MAP FILE FOR VARIOUS COMPOSITING OBJECTIVES 


When a cellular map file grows to a number of topics on resources 
dnd socio-economic factors, it readily lends itself to various 
compositing objectives. Typical objectives are shown to the right 
Of the map file, which may be composited with logical weights. 

This quick referencing is useful in environmenta.1 impact studies, 
growth area simulation for utilities, transportation, recreation 
etc. 
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subdivided into the dominant type cellular system and the 
percentile type cellular system. It is extremely important to 
ensure that the cell size used by the cellular system is able 
to provide adequate accuracy. 

Programs written on the cellular approach tend to be 
shorter and easier to write. This reduces the development costs 
of the software, as well as the computer cost of running each 
program. Versatility is also an advantage; the cellular approach 
allows the addition of new programs more easily and at ’less ex- 
pense. Some of the most manipulative facilities of thb' cellular 
format are: 

e Conversion of any map legend or data scale into other 
symbols, reranking, ordinal ranking, etc. without 
altering the original stored symbol. 

9 Computation of frequencies of the original symbols in 
the total map areas or in the subareas, via a master 
mapping procedure for setting up the subunits for 
automated mapping of any data so cataloged. 

« Logical selection to be applied to any pair of maps 
and carried to the next, etc., essentially using one 
map as a constraint upon the next, 

e More complex relationships betxveen a "dependent map" 
and a number of suspected independent maps may be dis- 
closed through statistical programs of multiple regression, 
discriminant analysis or any chosen program for regional 
phenomena. This may answer an analytical problem or 
may supply relationships for a new simulation composite 
map. 

In this region, the following applications of cellular map 
compositing have been made (see Appendix B for further information 
on applications) : 

® Urban development feasibility composites. 

• Stripmining reclamation feasibility. 

» Composite mapping revegetation scenarios. 

9 Power plant siting feasibility. 

• Selecting parkland in a metropolitan county. 

• Defining areas of housing restrictions in a floodable 
metropolitan county. 

• Anticipating critical areas of land use conflict over 
an entire state. 

• Conducting industrial optimum location search over a 
four-state region. 

• Simulating changes in riverine vegetation and water 
conditions under various management alternatives. 

• Allocating hospitals and field health services in a 
state according to composite indicators of local need. 
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6.1.2 Multisource Data 


Most spatial information, in whatever original form, 
can be converted into cellular form, as illustrated in the 
following Figure 6.1.2. Some of the more common grid conversions 
may be described as follows: 

• Any conventional map (such as topography, culture, 

or thematic patterns) may be optically projected upon 
a grid of appropriate fineness, encoded by one of 
several processes, and turned into data arrays 
stored in the cellular map file. A useful, visible 
beginning procedure for starting a cellular map file 
is to use manual mapping, encoding and keypunching; 
but larger jobs may use polygonal encoding or direct 
drum (raster) scanning. For example, air photo maps, 
after rectification and thematic photo interpretation, 
may be converted into grids by any of these procedures. 

• Another common form of input is data tables or tapes 
for preset areas, census tracts, zip code zones, land 
ownership jurisdiction, etc. for which a master map 
with given area codes is created in cellular form, 
after which tabular data values are automatically 
located via these location codes. 

• Another form of input is random sampling such as water 
well data, forest condition data, or climatic data. 

These data points, falling in certain cells, are inter- 
polated into a smooth surface by certain subroutines 

in cellular mapping. 

• With the general increase in available hardware for 
digitizing polygonal patterns, spatial engineering 

data and aerial photo interpretation are being digitized. 
The original polygonized format for national land use 
will become operational through the USGS LUDA mapping 
program. Then, in order to utilize various polygonal 
inputs for the purposes of uniform cellular mapping and 
analytical compositing among numerous maps in a data bank, 
the polygonal tape records can be read into the cellular 
format . 

• In general, for regional data banking, using the LANDSAT 
picture elements or larger grid cells and a series of 
other topics, the grid format will be versatile and low 
in cost. This is further illustrated by Figure 6.1.3. 
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FIGURE 6.1.2 
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FIGURE 6.1-3 
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- The valuable role of the LANDSAT CCT inputs in a cellular 
map file is to provide the current land use and cover information, 
making ivorthwhile the assembly of numerous other data on resources 
and socio-economic factors, and offering rapid compositing service 
to various state and federal agencies, counties, etc, 

6,lo3 Project Software 


In developing this particular project, several versions of 
cellular mapping and procedure in compositing were advanced. An 
earlier version of a computer composite mapping system (CMS-Ij 
was introduced into the region by the U.S. Department of Commerce 
via a 1970 economic study by four states: Colorado, Utah, New 

Mexico and Arizona, This enabled them to make an optimal industry 
location search over a large region for some 30 industries. Due to 
its peculiar language configuration which limited the program to cer- 
tain computers, the original program was reprogrammed and features 
were added by the Federation — this became the CMS-II program. In 
the meantime, two universities developed special versions of their 
o^m: (1) the Resources Analysis Program (RAP) of the Bureau of 

Business Research at the University of Utah, which uses polygonal 
inputs and a Calma digitizer; and (2) the Generalized Map Analysis Pla: 
ning System (GMAPS) program of the Colorado School of Mines which 
runs in an interactive mode for convenient student utilization 
and fits onto a smaller core memory machine. 

The GMAPS version of cellular mapping was utilized in this 
project by the Los Alamos Scientific Laboratories for compositing 
the various maps prepared by four states — Arizona, Colorado, Ne;>r 
Mexico, and Wyoming (see Appendix B) . Montana and Utah then utilized 
their own systems for compositing (see Appendices E and G) . 

In the Energy Planning Division of the State of Montana, 
still another version of cellular mapping, the Environmental Re- 
sources Geo- Information System (ERGIS), l^ras developed to run as 
part of a drum raster scanning process. 

Thus, a variety of new applications of cellular mapping and 
integration of LANDSAT CCT maps were initiated in the region. Fig- 
ure 6.1.3. shows the conceptual linkage of LANDSAT and multisource 
data files. 
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6 . 2 Procedures 


6,2.1 LASL Compositing 


The states of Arizona, Colorado, New Mexico, and 
Wyoming supplied data for computer compositing to LASL. 

Each state collected the desired ancillary data, coded 
that data onto a 1.1 acre representation cellular grid (see 
Figure 6.2.1), and punched that information onto computer 
cards. CSU supplied the already digitized LANDSAT data. 

This "package" was then sent to LASL for processing. 
Supplier-processer interactions included editing and 
assigning values or weights to the data. 

* This project utilized for compositing the Generalized 
Mapping Analysis Planning System (GMAPS) implemented at Los 
Alamos in 1975. The system is functionally identical to 
the GMAPS program developed at the Colorado School of Mines. ^ 

I 

Computer-aided versions of composite mapping represent 
an important extension of traditional manual graphic procedures. 
The use of a computer requires some modification to the pro- 
cedures so that the mapped source data can be entered into 
computer storage. This process is called "map digitization." 
Source information is quickly digitized if the information is 
first plotted on a base map of suitable scale; then the map 
is overlaid by a coordinate grid and the conditions within 
each grid cell recorded. Not all cells need to be filled in 
by hand. The GMAPS program . for a sequence of similar 
cells to be "filled in" by the computer which scans the lines 
from left to right. The computer can also duplicate identical 
lines. On an average, only about one- third of the cells 
needs to be hand coded. The GMAPS program performs a series 
of logical checks and issues appropriate error messages. 

In general, this procedure is most efficient if arbi- 
trary character codes are assigned to each condition portrayed 
on the map. No connotation of relative value is intended at 
this point. For example, areas coded "A" are not necessarily 
better or worse for any particular use than areas coded "B." 
Subsequently, as shown in Figure 6.2.1, these matrix- analog codes 
can be converted to numerical values reflecting the capability 
of a particular coded condition' to support a particular activity. 
Such valuation procedures are performed by the computer very 
rapidly, allowing a single data component to have a number of 
different patterns of values, as appropriate to the analysis. 


^Turner, A.K. , "Computer Aided Environmental Impact Analysis Parts 
1 and 2," Mineral Industries Bulletin , Vol, 19, Nos. 2 and 3, 1976, 
Colorado School of Mines. 
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8. RESULTANT LAND PATTERN 

Site capacities from analysis maps can 
be easily visualized using a graaation of 
tones from dark to light. 


7. ANALYSIS MAP 

*As the data Is combined, descriptions 
of the site's capacity can be 
determined, -mapped, and quantified. 
The diagram-shows a hypothetical case 
of timber productivity with "6" 
showing high timber productivity and 
showing low productivity. 

a OTHER DATA MAPS 

Information on many resource variables 
IS known for each grid cell and each is 
stored as a separate data file. 

5. COMPUTER DATA DISPLAY MAP 

Data can be displayed on computer 
printed maps as symbols or grey tones. 


4. COMPUTER CARDS 

Coded data is keypunched onto compu- 
- Ter cards. Positions on computer cards 
correspond to those on the coding 
sheets. 

3. DATA CODING 

The coding process involves the 
determination of the predominant 
resource class within each grid ceil and 
recording it with a numerical symbol. 
The example in the diagram shows: 
1 = grassland; 2 = marsh, and 3 = forest. 

2. GRID 

The scale of the base map and the size 
of the grid cell are established 
according to the detail of the data 
available, and the detail of desired 
results. Thus, cell size could range from 
one acre to one square mile. 


1. SITE 

A site Is composed of many natural 
land features including soil, water, and 
vegetation. Each of these is identified 
3nd record<>d on j mao. 
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The GMA.PS program allows the overlaying, or "compositing” of 
maps, using arithmetic and logical expressions. "Arithmetic compositing" 
is a simple extension of the tonal overlay procedures; but it allows 
for much more varied analyses using combinations of additipii, sub-- 
traction, multiplication, and -div-i-s-ion in conjiinction with the ability 
to weight some components more heavily than others. 

The GMAPS prograni differs from the earlier programs ..be'cause 
it is designed to interact from a teletype terminal in a time-sharing 
environment with a computer. This makes GMAPS very attractive to use 
because: 

• the program asks a sequence of questions to which the user 
responds,' thereby defining the operations and sequence of 
operations he wishes the computer to perform; 

• the program allows the user to verify and correct commands, 
so that meaningless operations are eliminated; 

• the system is easily used by laymen; 

• the time-sharing concept gives the user access to a high- 
capacity computer at economical cost. 

In making modifications to the original GMAPS program, LASL 
specified continued corapatability with the smaller machine (DEC-10) 
versions. This proved to be a very real advantage, not merely a theoreti- 
cal one. For this project, initial compositing runs were begun at Los 
Alamos, but many of the composites were completed at the Colorado School 
of Mines. 

The LASL versions of GMAPS are being ejqjanded to take advantage of 
the larger computing poiver and more sophisticated data entry_and display 
hardware available at the Laboratory. For this project, the ability to 
directly generate 35mm colorslides of the maps and composites was demon- 
strated to the state agencies. IVhile LASL is not in the position to 
offer such services on a routine basis, several commercial firms can do 
so, and it therefore seemed appropriate to demonstrate these capabilities. 

Interfacing procedures were required to convert the LANDSAT data 
received from CSU to a form compatible with GMAPS. Figure 6.2.2 shows how 
the interfacing was accomplished. There were, for purposes of conven- 
ience, two computer programs called ERTSl and ERTS2. 

'The first program, ERTSl, accepted the CSU tapes and through 
interactive questioning to the user converted the five character CSU 
land-use codes found in the first tape record into easily understand- 
able legend information. The remaining data on the CSU tape were 
divided into two files containing formated record sequences. These 
records were designed to allow them to be data cards. One file 
sequence contained the legend information and the symbol strings 
corresponding to the land-use definitions for each cell (or pixel) . 

The other file contained the discriminant values for each cell. Both 
files contained one record set for each line across the entire map 
width. In other words, each file contained a matrix in row-ordered 
form. 
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The second program, ERTS2, converted these two intermediate 
. files into a single GMA.PS input file. These files are organized into 
"sectors," a sector being a block of 120x120 cells.. Within each 
sector, two sequences of 60-character records, corresponding to the left 
and right halves of the sector, are placed in sequence. These file 
structures are termed the "P-cards" format and form the common data 
source formats for GMA.PS and CMS program types. 

ERTS2 performed one other important function, the thresholding of 
the data. At the user's option, points with a low statistical signifi- 
cance may be thresholded out by comparing their discriminant value to 
the Chi-square value corresponding to the nxiraber of channels used in 
the classification and the desired threshold level. The procedure is 
identical to that defined in the RECOG system part of LMS) documenta- 
tion.^ If the discriminant value exceeded the Chi-square value, ERTS2 
substituted a predefined "threshold" character for the existing classifi- 
cation. 


This process allows for the generation of a large number of 
different thresholded variations of the basic LANDSAT derived land- use 
map data. This gives maps of various accuracies for certain land-use 
types. Each such product can be handled by GMAPS as a separate input 
source. It was primarily for this reason that the interfacing opera- 
tion was divided into t\ro stages, corresponding to ERTSl and ERTS2. 

It should be noted that both steps in the process were designed 
as interactive procedures, allowing direct interface \^ith the user via 
the teletype terminal. 


6.2.2 Montana Compositing 

The Energy Planning Division (EDP) of Montana's Department of 
Natural Resources and Conservation as lead agency conducted the ancillary 
data collection, digitizing and compositing. CSU again supplied the 
LANDSAT digitized data, 

EDP's equipment includes a raster scanner which reads "data" 
directly from maps by shining a light on the document and measuring the 
intensity of the reflected light. These light readings vary by color 
and therefore provide a "code" similar to the cellular process. 

The new data is written to a magnetic tape for computer process- 
ing. 


The created map file had the same cell size and registration as 
the LANDSAT data. Compositing of the multisource map files was then 
conducted . 


^T. Ells, L.D, Miller, and J.A. Smith, User's Manual for RECOG CPattern 
RECOGnition Programs ), Science Series No. 38, Department of Watershed 
Science, College of Forestry and Natural Resources, Colorado State Univer- 
sity, Fort Collins, Colorado. 
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6*2.3 Utah Compositing 


Utah used a polygonal approach on a Calma digitizer, tape 
files therefrom to be applied in a new Resource Analysis Procedure (RAP) 
program developed by the Bureau of Economic and Business Research at the 
University of Utah. This is a second generation beyond CMS, allowing in- 
put either as cells or polygons. 

Utah utilized the RAP to process the LANDSAT imagery vin conjunction 
with other factor maps pertinent to the Farmington quad. The RAP program 
has gone beyond the CMS based systems in that it will permit input 
from various other media, and then internally make them all compatible. 

The reason for this was to broaden the input data bases for factor maps 
and eliminate some of the restrictions which are imposed by a strictly 
cellular ized, card-oriented system. For example, by implementing a 
polygonal approach in RAP, steps 2, 3 and 4 as shown on figure 6.2.1 
may be eliminated and niaps may be produced at any scale desired. 

Polygonal input processing routines within RAP can handle maps digitized 
from several different projections (i.e., Lambert, U.T.M., etc.) and 
internally translate them to geographic coordinates for compatibility. 
Since RAP can maintain maps in a real-valued state, mathematical manipu- 
lation for compositing purposes has true meaning, whereas it does not 
for maps having ordinal, numerical values. A real-valued map would 
contain the actual value at each cell across a map topic, thereby per- 
mitting true compositing of characteristics and allowing for the genera- 
tion of a large number of different threshold variations. 

RAP’s compositing features enable the user to produce data maps 
as a function of many other topic maps. A new topic map (composite) 
can be produced through the use of logical expressions, recoding, 
thresholding, arithmetic expressions, boolean expressions, mathematical 
functions, and weighting. RAP has the feature of including entire maps 
or only component parts of various maps in this compositing process • 

These new composited topic maps can then again be functional variables 
in other compositions. Real-valued, ordinal and non-numeric map topics 
may be used in the same compositing process. 

This project gave the planners in Utah the opportunity to imple- 
ment the interface between CSU-produced LMS tapes and the RAP system. 

Along with this interface RAP will also interface with U.S. Forest 
Service COMLUP system files, USGS GIRAS (LUDA) polygonal data files, 

Bendix Corporation-produced LANDSAT imagery, IDS digitized data bases, 
SYMAP map files and the existing methods of the original CMS system. 

For the user’s convenience, RAP can be run in either batch mode 
or demand (interactive) mode. The latter operational mode enables the 
terminal operator to make real-time decisions concerning the manipula- 
tions of map topics. Demand mode is enhanced by immediate user awareness 
of error conditions so that meaningless commands can be corrected or 
eliminated. RAP will prompt the operator with appropriate messages 
when input is desired or various options present themselves. 
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RAP has maintained all the desirable features of the CMS system, 
such as aggregate area dictionary s3mibol conversion, map sector input 
via P-cards, internal and external symbol conversion, external data file 
input via master maps, boundary overlays, reduction of maps, shifting 
of origins, map interpolation with limited known points and statistical 
analysis* In addition a linear programming routine has been incorporated 
with the options of inputting data from topic maps for the purpose of 
optimizing spacial values. This feature, however, was not used in this 
-demonstration. 
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FIGURE 6.3.2 


ARIZONA COMPOSITE ANALYSIS 



Negating factors include flood hazard, nonprivate land, prime agricultural 
land, and severe soil limitations for residential construction. 

Current residential developments (from LANDSAT Topic Map) were classified 
as potentially available land, thereby allov^ing the suitability of such 
areas to be evaluated according to the site determining parameters. 
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6 . 3 Results 


Each state developed demonstration composi-ting- analyses for their 
comprehensive test sites, as shown below. 


FIGURE 6.3.1 


COMPREHENSIVE TEST SITE ANALYSIS 


State 

Comprehensive Test Site 

Computer Composite Analysis 

Arizona 

Hedgepeth Hills 

Optimal locations for resi- 
dential development 

Colorado 

Fox Creek 

Index of revegetation 
capabilities 

Montana 

Colstrip, SE 

Mined land reclamation 
feasibility 

New Mexico 

Santa Fe 

Restrictions against urban 
development 

Utah 

Farmington 

Constraints to nonagricul- 
tural development 

Wyoming 

Buffalo 

Physical limitations to 
urban growth 


In the following Figures, each state's compositing procedure is 
charted, followed by a reduced version of their final composite map. 

In the charts, input topics are indicated in the upper boxes, relative 
weights (values) below and the resulting composite in the lower box. 
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FIGURE 6.3.3 
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FIGURE 6.3.4 

COLORADO COMPOSITE ANALYSIS 
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FIGURE 6.3.5 
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FIGURE 6.3.6 


MONTANA COMPOSITE ANALYSIS 





Coal Mining 
Suitability 
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Coal Mining Suitability Composite Map (Figure 6.3.7, following page) 
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FIGURE 6 . 3.7 

Coal Mining Suitability Composite Map 
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Figure 16 Coal Mining Suitability Composite Map 
(Legend on next page) 
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FIGURE 6.3.8 


NEW MEXICO COMPOSITE ANALYSIS 



{ ) Map Number Reference In Appendix 3 

O Weights Assigned To The Data Source During Compilation 
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FIGURE 6.3.9 
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FIGURE 6.3.10 
UTAH COMPOSITE A^^ALYSIS 
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UTAH COMPOSITE MAP 


FIGURE 6.3.11 
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FIGURE 6.3.12 




WYOMING COMPOSITE ANALYSIS 



Only a portion of the complete test site was analyzed due to a lack of adequate 
data. 
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FIGURE 6.3.13 
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